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FOREWORD 


This  annual  report  for  the  continued  development 
of  a  dual-chamber  controllable  solid  propellant  rocket 
motor  (DCCSR)  describes  the  progress  during  the  third 
year  of  this  program,  which  is  sponsored  by  the  Air 
Force  Rocket  Propulsion  Laboratory,  Edwards  Air 
Force  Base,  California.  The  research  and  development 
efforts  of  the  program  are  being  performed  by  Northrop 
Carolina,  Inc. ,  a  Subsidiary  of  Northrop  Corporation, 
Asheville,  North  Carolina,  Tinder  Air  Force  Contract 
AF  04{61 1)-9067.  This  report  is  presented  in  two 
volumes:  Volume  I  -  Research  and  Development  Efforts, 
and  Volume  II  -  Analytical  Study.  This  volume 
(Volume  II)  presents  the  results  of  a  study  to  determine 
the  effect  of  motor  performance  parameters  and  propel¬ 
lant  characteristics  on  the  mass  fraction,  burnout 
velocity,  and  motor  envelope  of  the  dual-chamber 
controllable  solid  propellant  rocket  motor. 


NOTICE 


Northrop  Carolina,  Inc. ,  has  been  assigned  a  patent  application  by  the  U.  S. 
Patent  Office  to  cover  the  Controllable  Solid  Propellant  Rocket  Motor  invention 
disclosed  in  this  publication,  and  the  Commissioner  of  Patents  has  issued  a  secrecy 
order  thereon.  This  secrecy  order  requires  that  those  who  receive  a  disclosure  of 
the  subject  matter  be  informed  of  the  existence  of  the  secrecy  order  and  of  the 
penalties  for  the  violation  thereof. 

The  recipient  of  this  report  is  accordingly  advised  that  this  publication 
includes  information  which  is  now  under  a  secrecy  order.  It  is  requested  that  he 
notify  all  persons  who  will  have  access  to  this  material  of  the  secrecy  order. 

Each  secrecy  order  provides  that  any  person  who  has  received  a  disclosure 
of  the  subject  matter  covered  by  the  secrecy  order  is 

"in  nowise  to  publish  or  disclose  the  invention  or  any  material 
information  with  respect  thereto,  including  hitherto  unpublished  details 
of  the  subject  matter  of  said  application,  in  any  way. to  any  persons  not 
cognizant  of  the  invention  prior  to  the  date  of  the  order,  including  any 
employee  of  the  principals,  but  to  keep  the  same  secret  except  by  written 
permission  first  obtained  of  the  Commissioner  of  Patents,  under  the 
penalties  of  35  U.  S.  C.  (1952)  182,  186" 

Although  the  original  secrecy  order  forbids  disclosure  of  the  material  to 
persons  not  cognizant  of  the  invention  prior  to  the  date  of  the  order,  a  supplemental 
permit  attached  to  each  order  does  permit  such  disclosures  to: 

"(a)  Any  officer  or  employee  of  any  department,  independent  agency  or 
bureau  of  the  Government  of  the  United  States. 

(b)  Any  person  designated  specifically  by  the  head  of  any  department, 

independent  agency  or  bureau  of  the  Government  of  the  United  States, 
or  by  his  duly' authorized  subordinate,  „as  a  proper  individual  to  receive 
the  disclosure  of  the  above  indicated  application. 

The  principals  under  the  secrecy  are  further  authorized  to  disclose 
the  subject  matter  of  this  application  to  the  minimum  necessary  number  of 
persons  of  known  loyalty  and  discretion,  employed  by  or  working  with  the 
principals  or  their  licensees  and  whose  duties  involve  cooperation  in  the 
development,  manufacture  or  use  of  the  subject  matter  by  or  for  the  Govern¬ 
ment  of  the  United’  States,  provided  such  persons  are  advised  of  the  issuance 
of  the  secrecy  order.  " 

No  other  disclosures  are  authorized,  without  written  permission  from  the 
Commissioner  of  Patents.  Penalties  for  violation  of  a  secrecy  order,  include  a 
fine  of  up  to  $10,  000  or  imprisonment  for  not  more  than  two  years  or  both. 

It  must  be  understood  that  the  requirements  of  the  secrecy  order  of  the 
Commissioner  of  Patents  are  in  addition  to  the  usual  security  regulations  which 
are  in  force  with  respect  to  classified  materials  in  this  report.  The  usual  security 
regulations  must  still  be  observed  notwithstanding  anything  set  forth  in  the  secrecy 
order.  In  the  event  that  this  report  shall  be  declassified,  the  secrecy  order  remains 
in  full  force  until  it  is  specifically  rescinded. 
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SECTION  I  -  INTRODUCTION 


I.  GENERAL 

Northrop  Carolina,  Inc.  ,  a  subsidiary  of  Northrop  Corporation, 
has  been  developing  a  dual- chamber  controllable  solid  propellant 
rocket  motor  (DCCSR)  under  sponsorship  of  the  Air  Force  Systems 
Command,  Rocket  Propulsion  Laboratory,  Edwards  Air  Force  Base, 
California.  This  program,  now  in  its  third  year,  is  funded  under 
Contract  AF  04(61 1)-9067.  The  progress  of  the  program  has  been 
reported  in  quarterly  and  annual  technical  documentary  reports 
(References  1  through  13).  This  volume  presents  the  results  of  a 
separate  parametric  study  conducted  as  part  of  the  third  year's 
effort. 

The  DCCSR  concept  is  fully  described  and  illustrated  in 
Volume  I.  However,  a  brief  description  is  given  here  to  provide 
a  better  understanding  of  the  information  and  data  presented 
in  this  volume.  The  concept  utilizes  two  propellant  chambers  separated 
by  a  throttle  valve.  The  forward  chamber  contains  a  cool-burning 
fuel-rich  propellant;  the  aft,  an  oxidizer-rich  propellant.  The  fact 
that  the  forward- chamber  combustion  products  are  relatively  cool 
burning  (2800°F)  permits  state-of-the-art  materials  to  be  used  on 
the  aft  chamber,  where  high-temperature  combustion  takes  place. 

A  multiple  pyrogen  ignition  system  is  included  as  a  part  of  the  forward 
chamber.  (Confidential) 

Thrust  is  initiated  by  igniting  the  forward  propellant  with  a 
single  pyrogen.  The  relatively  cool  combustion  products  from  the 
forward  chamber  are  thrbttled  through  the  control  valve  into  the  aft 
chamber  where  additional  thermochemical  reaction  occurs,  resulting 
in  more  energy  release.  The  aft  propellant  will  not  burn  without  the 
heat  supplied  by  the  forward  propellant.  Throttleability  is  achieved 
by  varying  the  forward- chamber  pressure  (and  burning  rate)  by  varying 
the  position  of  the  valve.  The  aft  propellant  actively  supports  com¬ 
bustion  when  the  combustion  gases  from  the  forward  propellant  pass 
over  it.  (Confidential) 
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Thrust  can  be  terminated  at  any  time  during  the  burning  period  by 
suddenly  increasing  the -valve  flow  area,  which  produces  a  rarefaction 
wave  that  extinguishes  combustion  of  the  forward  propellant.  Since  the 
aft  propellant  will  not  sustain  combustion  except  at  high  chamber 
pressures  and/or  temperatures  above  300  C  without  an  external  heat 
source,  it  too  is  extinguished.  The  on-off  cycle  operation  can  be 
repeated  on  command  by  reigniting  the  forward  propellant  using  another 
pyrogen  igniter  of  the  multiple  ignition  system  for  each  restart.  (Confidential) 

2.  PURPOSE  OF  STUDY 

This  study  was  conducted  in  order  to  establish  trends  and  optimum 
conditions  of  mass  fraction,  boost  velocity,  and  motor  envelope  with 
variations  in  motor  performance  parameters  and  propellant  character¬ 
istics  for  the  dual-chamber  controllable  solid  propellant  rocket  motor 
concept.  It  should  be  noted  that  the  results  presented  herein  do  not 
necessarily  represent  optimum  designs;  rather,  off-optimum  as  well  as 
optimum  conditions  have  been  examined. 

3.  SCOPE 

Five  independent  variables  (motor  performance  parameters  and 
propellant  characteristics )were  selected  for  this  study:  (1)  total  impulse, 

(2)  minimum  thrust,  (3)  thrust  throttling  range,  (4)  number  of  on-off 
cycles  available,  and  (5)  motor  specific  impulse.  The  first  four 
variables  define  the  capability  of  a  throttleable  stop-.restart  motor,  and 
each  directly  affects  mass  fraction,  boost  velocity,  and  motor  envelope. 

The  fifth  variable,  which  is  a  function  of  propellant  composition, 
directly  affects  motor  size  and  boost  velocity. 

The  range  of  total  impulse,  minimum  thrust,  and  thrust  throttling 
range  values  included  in  this  study  are  presented  in  Table  I.  Minimum 
thrust  was  varied  from  l/500  to  l/20  of  total  impulse,  that  is,  for 
maximum  burn  times  of  20  to  500  sec.  The  thrust  throttling  ranges 

considered  were  1  to  1, _ 5  to  1,  and  20  to  1,  and  the  number  of  on-off 

cycles  used  was  1,  10,  20,  and  40.  Vacuum  specific  impulse,  at  an 
expansion  ratio  of  20  to  1,  was  varied  from  265  to  280  to  300  lb-sec/lb. 

It  was  felt  that  these  ranges  of  independent  variables  would  encompass 
those  required  in  most  applications  of  a  DCCSR.  (Confidential) 

The  analysis  was  divided  into  two  phases.  In  Phase  I,  the 
independent  variables  listed  above  were  investigated  separately,  with 
propellant  properties  and  structural  materials  fixed.  In  Phase  II,  case 
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TABLE  I  -  TOTAL  IMPULSE  AND  THRUST  RANGES 


FOR  PARAMETRIC  ANALYSIS 


Maximum  Thrust  (lb.) 


Case  1 


200 
500 
1,  000 
5,  000 

1,  000 
2,  000 
5,  000 
25,  000 

2,  000 
4,  000 
10, 000 
50,  000 


Case  2 


100 
250 
500 
2,  500 

1,  000 
2,  500 
5,  000 
25,  000 

5,  000 
10, 000 
25,  000 
125,  000 

10*000 
20,  000 
50,  000 
250, 000 


Case  3 


Mil 

rm 


10, 000 


10, 000 

2 


100, 000 


40, 000 
100, 000 
500, 000 

40, 000 
80, 000 
200, 000 
1,  000,  000 


(Confidential) 
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material,  propellant  mixture  ratio,  burning-rate  constants,  and 
burning-rate  pressure  exponents  were  varied  as  follows: 

1>  The  effect  of  three  case  materials  (steel,  titanium  and 
fiberglass)  was  evaluated  for  (1)  the  four  thrust  levels 
in  the  100,  000-lb-sec  motor,  (2)  the  other  impulse 
levels  at  a  minimum  thrust  corresponding  to  200-sec 
operating  time,  and  (3)  for  throttling  ratios  of  1,  5, 
and  20. 


2.  The  effect  of  aft-to-forward-chamber  propellant  mixture 
ratios  (2  to  1,  3  to  1,  and  4  to  1)  was  investigated  over 
the  total  impulse  and  minim um  thrust  range  at  a  constant 
throttling  ratio  of  5  to  1. 


3.  The  effect  of  burning- rate  constants  was  investigated  for 
the  four  thrust  levels  in  the  100,  000-lb-sec  motor  at  a 
throttling  ratio  of  5  to  1. 


4.  The  effect  of  propellant  burning-rate  (pressure)  exponent 

was  evaluated  at  throttling  ratios  of  1  to  1,  5  to  1,  and 

20  to  1  at  a  minimum  thrust  of  500  lb^  in  the  100,  000-lb-sec 

motor.  The  burning-rate  exponents  evaluated  were: 

(1)  for  the  forward  chamber,  0.6,  0.8,  and  0.9;  (2)  for 

the  aft  chamber,  0.8,  1.0,  and  1.1.  ...  ..  ,, 

(Confidential) 

The  range  of  variables  and  constants  used  in  this  study  are 
listed  in  detail  in  Section  II,  paragraph  1. 


4.  METHOD  OF  COMPUTATION 

This  study  was  carried  out  by  means  of  Northrop  Carolina's  IBM 
1620  Data  Processing  System  and  an  IBM  1622  Card  Read  Punch,  which 
assured  accurate  and  rapid  processing  of  the  required  data.  The  data 
processing  system's  limited  storage  capacity  (20,000  digits)  neces¬ 
sitated  the  processing  of  data  through  four  separate  computer  sub¬ 
routines.  These  subroutines,  prepared  by  Northrop  Carolina,  con¬ 
sisted  of  (1)  a  steady-state  internal  ballistics  subroutine,  (2)  a  grain 
configuration  subroutine,  and  (3)  two  weight  subroutines.  The  infor¬ 
mation  flow  through  these  subroutines  is  shown  in  Figure  1.  A  detailed 
description  of  the  subroutines  and  the  overall  computer  program  is 
presented  in  Section  II. 
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SECTION  II  -  TECHNICAL  APPROACH 


1.  PARAMETER  SELECTION 
a.  General 

As  mentioned  in  Section  I5  this  parametric  study  was 
conducted  in  two  phases.  The  effect  of  the  major  ballistic 
parameters  was  investigated  in  Phase  I;  the  remaining,  less 
significant  parameters  in  Phase  II.  As  a  result,  more  com¬ 
puter  runs  were  required  for  Phase  I  than  for  Phase  II. 

Table  II  summarizes  the  parameters  investigated  in 
each  phase,  showing  which  were  varied  and  which  held  constant. 
Note  that,  as  shown  in  Table  II,  Phase  II  was  subdivided  into 
Phases  Ha  through  lid,  depending  on  the  parameters  that  were 
held  constant,  as  follows: 


Phase 

Variable 

II  a 

Case  material 

II  b 

Burning  rate 

II  c  Propellant  burning  rate  exponent 

II  d  O/F  ratio,  0  (theta) 

The  values  and  ranges  of  the  parameters  investigated  in 
each  phase  are  listed  in  paragraphs  b  and  c,  below,  respectively. 

b.  Phase  I  Parameters 

The  values  of  the  constants  used  in  Phase  I  are  listed  below. 


Constant  Value 

Case  material  Steel 

O/F  ratio,  0  3.  0 

Forward  propellant  density  (lb^/in.  3)  0.  053 

Aft  propellant  density  (lb^/in.  3)  0.  070 


AFRPL-TR-65-209,  Vo!  H  CONFIDENTIAL 

TABLE  II  -  PARAMETERS  INVESTIGATED  IN  EACH  PHASE* 


Parameter  — 

Phase 

I 

II  a 

II  b 

II  c 

II  d 

Total  impulse,  vacuum  (lb^-sec) 

V 

V 

C 

c 

V 

Specific  impulse,  vacuum  (Ib.-sec/lb  1 

t  m 

V 

C 

C 

c 

c 

Minimum  thrust  (lb^) 

V 

V 

c 

V 

Thrust  ratio,  maximum/minimum 

V 

c 

Starts 

V 

c 

Case  material 

C 

V 

c 

Aft  propellant  rate  constant  (in.  /sec) 

C 

C 

c 

Forward  propellant  rate  constant  (in.  /sec) 

C 

C 

c 

Aft  propellant  burning- rate  exponent 

C 

C 

C 

c 

Forward  propellant  burning- rate  exponent 

C 

C 

C 

c 

O/F  ratio,  0 

C 

c 

c 

c 

V 

Aft  propellant  density  (lb  /in.  ) 

m 

C 

c 

c 

c 

c 

Forward  propellant  density  (lb  /in.  3) 

C 

c 

c 

c 

c 

Aft  propellant  flame  temperature  (  F) 

C 

c 

c 

c 

c 

Forward  propellant  flame  temperature  (°F) 

C 

c 

c 

c 

c 

Aft  chamber C,  (sec  x) 

U 

c 

c 

c 

c 

■ 

Forward  chamber  C  ,  (sec  ) 

a 

C  i 

c 

c 

c 

c 

Minimum  aft  chamber  pressure  (psia) 

c 

c 

c 

c 

c 

Minimum  forward  chamber  pressure  (psia) 

c 

c 

c 

c 

c 

Expansion  ratio 

c 

c 

c 

c 

c 

I*  - ,  ■■ 

C  and  V  denote  constant  and  variable,  respectively. 

(Confidential) 
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Constant 

Forward  propellant  exponent,  n 

Aft  propellant  exponent,  n 

Forward  propellant  rate  constant  (in.  /sec) 

Aft  propellant  rate  constant  (in.  /sec) 

Forward  propellant  flame  temperature  (°F) 

Aft  propellant  flame  temperature  (°F) 

C^,  forward  chamber  (sec  l) 

C  ,,  aft  chamber  (sec  x) 
a 

Minimum  forward-chamber  pressure  (psia) 
Minimum  aft-chamber  pressure  (psia) 
Expansion  ratio 


Value 
0.  80 
1.  00 
0.  0008 
0.  0004 
2700 
5500 
0.  00765 
0.  00643 
100.  0 
50.  0 
20.  0 


(Confidential) 

The  range  of  the  variables  used  for  Phase  I  are  as 
follows: 


Variable 

Total  impulse  (lb^-sec) 
Maximum/minimum  thrust  ratio 
Minimum  thrust  (lb^) 

Specific  impulse,  vacuum  (lb^-sec/l* 
Starts 

Phase  II  Parameters 


Range 


104 

to  10s 

1. 

-WO 

,  000 

280,  300 

1,  10,  20 
(Confidential) 


The  less  significant  parameters  (case  material,  burn  rate, 
pressure  exponent,  and  theta)were  varied  in  Phase  II.  The  range 
of  these  variables  and  the  values  of  the  constants  used  in  each 
portion  of  Phase  II  are  listed  below. 


Phase  II  a 

Constants  -  Same  as  Phase  I,  except  for  the  following: 


-9- 


CONFIDENTIAL 


AFRPL-TR-65-209.  Vol  II C6NFI DENTI.il 


Case  Materia] 

Listed  below 

\ 

Specific  impulse,  vacuum 

_ 

[  '  I 

(lb  -sec/lb  ) 

X  rr. 

280 

Starts 

20 

L  -■ 

Variables 

Maximum/minimum  thrust  ratio 

1,  5,  20 

;  ’  I 

i 

L  /  * 

| 

Case  material 

Steel, 

i 

| 

Total  impulse  (Ib^-sec) 

Fiberglass, 

Titanium 

104  to  10® 

j 

f  i 

U 

Minimum  thrust  (lb^) 

50  to  5000 

Phase  11  b 

Constants  -  Same  as  Phase  1,  except  fo 

Burn  rate  constants 

r  the  following: 

Listed  below 

i  ! 

I.J 

;  i 

Specific  impulse,  vacuum 

(lb,- sec/lb  ) 
f  m 

280 

u 

II 

Starts 

20 

u 

« 

Maximum/minimum  thrust  ratio 

5 

t 

5 

1  . 

Total  impulse  (lb^-sec) 

106 

....  .Li 

Variables 

Aft  propellant  burn  rate 
constant  (in  /sec) 

0.  0001  to 

|  j 

L  J 

!  1 
f  > 

*  ? 

J 

• 

0.  0010 

Forward  propellant  burn  rate 
constant  (in.  /sec) 

0.  0002  to 

!  ( 

i  i 

0.  0032 

1 

Minimum  thrust 

1  200  to  500 

|  | 

Phase  II  c 

Constants  -  Same  as  Phase  I,  except  for  the  following: 

U 

1 1 

Burn  rate  exponents 

Listed  below 

LI 

Specific  impulse,  vacuum 

(lb,-  sec/lb  ) 
f  m 

280 

u 

Starts 
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2.  INTERNAL  BALLISTICS  SUBROUTINE 

The  internal  ballistics  subroutine  was  prepared  to  calculate 
ballistic  design  and  performance  parameters  for  the  DCCSR  motor 
under  steady- state  operating  conditions,  This  subroutine  contains 
two  principal  options  which  permit  the  program  to  be  used  to  design 
aft-chamber  parameters  and  calculate  performance  for  test  motors 
(Option  I)  or  to  design  motors  to  meet  specified  thrust  and  total 
impulse  requirements  (Option  II),  The  general  subroutine  logic  is 
shown  in  the  diagram  in  Figure  2;  the  actual  IBM  1620  computer  sub¬ 
routine  listing  is  given  in  Appendix  A. 

The  equations  employed  in  this  subroutine  are  conventional 
ballistic  relationships  which  have  been  modified  for  dual-chamber 
motors,  as  outlined  previously  (Reference  1).  The  forward 
chamber,  with  choked  flow  at  the  valve,  operates  as  a  conventional 
rocket  motor.  Aft-chamber  parameters  are  calculated  as  a  function 
of  the  ratio  of  aft-to-forward  mass  flow,  6  ,  which  under  steady- 
state  conditions  is  defined  as: 

(P  SaPn) 

e  - - sA.  ,  (1) 

<PSapn>fwd 

where 

p  -  propellant  density, 

S  -  burning  area, 
a  =  burn-rate  constant, 

P  ~  operating  pressure,  and 
n  =  pressure  exponent 


The  aft-chamber  surface  area  was  calculated  by  rearranging 
this  equation,  as  follows: 


Also  referred  to  as  burning- rate  exponent  herein. 
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Start 


Select  Option  and 
Read  Input 


Opticrrr  Number 


Compute  I  and  C_ 
forSp9  D 


Compute  Propellant  Weights,  Propellant 
Webs,  Propellant  Surface  Areas,  Aft 
Throat  Area,  and  Maximum  Burning  Time 


Compute  Aft  Chamber  Pressure 
for  Given  Forward  Chamber  Pressure 


Compute  Weight  of  Aft  Propellant, 
Forward  Throat  Area, 

Aft  Web,  and  Total  Impulse 


Inclement 

Forward 

Pressure 


Go  to  Next 
Subroutine 


Figure  2  -  Internal  Ballistics  Subroutine  J-iOgic  mag  ram 


r 
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The  over-all  mass  flow  balance  for  the  motor  is 


M,  ,  +  M  .  =  M  =  (C  ,  P  A  )  ,  , 
fwd  aft  tot  '  d  t'aft 


(3) 


where 


Cd  =  discharge  coefficient  at  aft  nozzle,  and 


A^  =  aft  nozzle  throat  area. 


Hence, 


or 


<PSaAwd  +  <  PS^Kit=  <CdPAt  >aff 


1  +  0  = 


^cdPAtU 

^SapI1>fwd 


(4) 

(5) 


The  aft-chamber  throat  area  was  calculated  by  rearranging 
Equation  5;  that  is, 


n 


(pSaP  .  ).  , 

, ,  .  fll  min  fwd 

\ft  =  <1+9>  «=dPmJaf, 


(6) 


Aft-chamber  pressure  was  calculated  as  a  function  of  forward- 
chamber  pressure  from  Equation  4  by  an  iterative  procedure.  Initially, 
an  estimate  of  aft-chamber  pressure  was  substituted  into  the  left  side 
of  Equation  4,  and  the  aft-chamber  pressure  on  the  right  side  was 
calcui.-iod  from 


..  ^SapnLd  * 


aft 


<cd  \> 


aft 


'rl>  cl'  .culated  pressure  was  then  compared  with  the  estimated 
piujsare.  If  the  two  did  not  agree  within  0.  5  percent,  another  estimated 
pxessure  was  obtained  by 

^est  ^  ^calc.  ^est, 

n  (n-1)  (n-1) 

and  the  procedure  repeated  until  convergence  was  achieved. 
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When  the  calculated  pressure  agreed  with  the  estimated  pressure 
within  0.  5  percent,  0  was  calculated  frore  Equation  1  and  compared 
with  the  estimated  0  from  which  C  ,  was  obtained.  If  these  did  not 

daft 

agree  within  0.  5  percent,  a  new  0  was  estimated  and  aft-chamber 
pressure  was  calculated  using  a  new  C  .  .  This  procedure  was 

aft 

repeated  until  convergence  was  achieved.  Thrust  was  then  calculated 
from  the  product  of  specific  impulse  at  the  final  0  value  and  mass 
flow  rate. 

3.  GRAIN  CONFIGURATION  SUBROUTINE 

The  grain  configuration  subroutine  was  prepared  for  determining 
the  most  suitable  grain  configurations  for  the  DCCSR  forward  and 
aft  chambers.  The  subroutine  calculated  the  diameter,  length,  port 
area,  and  volumetric  loading  of  the  grains,  but  did  not  calculate 
the  detailed  grain  design.  The  four  basic  grain  types  employed,  in 
order  of  increasing  mass  flow  requirements,  were  end  burning,  cylin- 
drically  (center)  perforated,  starr  and' wagon-wheel'  designs.  "  The S e 
configurations  were  selected  for  either  a  cylindrical  or  spherical  forward 
chamber  ar*d  for  a  cylindrical  aft  chamber,  .  either  in  tandem  or  around  the 
outside  of  the  forward  chamber.  The  required  input  for  this  sub¬ 
routine,  supplied  by  the  steady-state  internal  ballistics  subroutine, 
consisted  of  forward-  and  aft-grain  burning  surface  areas  and  web 
thicknesses,  ,  and  aft-chamber  throat  area.  The  subroutine  output 
provided  the  dimensions-neceasary  for  case  and  insulation  design. 

The  subroutine  logic  is  shown  in  the  diagram  in  Figure  3,  while  the 
actual  computer  subroutine  listing  is  given  in  Appendix  A. 

The  equations  employed  in  this  subroutine  are  general  geometry 
equations  that  relate  diameters  to  areas  and  volumes  of  both  cylinders 
and  spheres;  the  actual  equations  used  are  shown  in  the  program 
listing  for  this  subroutine  in  Appendix  A.  The  subroutine  initially 
established  the  forward-grain  configuration  and  its  dimensions,  based 
on  forward-grain  surface  area  and  web  thickness  inputs.  First,  a 
cylindrical  end-burner  design  was  assumed  and  the  grain  diameter 
calculated  from  the  inputs.  If  the  length-to- diameter  ratio  of  this 
design  was  less  than  0.  5,  the  cylindrical  end-burner  design  was 
rejected  and  a  spherical  end  burner  assumed.  The  program  thus 
proceeded  from  designs  with  high  web  fractions  to  those  with  low  web 
fractions  until  a  design  with  a  reasonable  length-to -diameter  ratio 
was  obtained.  The  program  then  printed  out  the  type  of  configuration, 
grain  outside  diameter,  length,  and  volumetric  loading.  .  For  cylindrical 
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Start 


Figure  3  -  Grain  Configuration  Subroutine  Logic  Diagram 
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center  perforated  and  star  designs,  the  program  would,  when  feasible, 
give  an  alternate  spherical  forward-chamber  design. 

After  the  forward  grain  configuration  was  selected,  the  subroutine 
had  the  option  of  proceeding  to  another  case  or  establishing  the  aft- 
grain  configuration.  The  aft-grain  program  was  similar  to  the  forward 
one  in  that  designs  with  increasing  web-to1- diameter  ratios  were 
considered  until  a  design  giving  a  reasonable  length-to-  diameter  ratio 
was  achieved.  An  attempt  was  made  to  design  an  aft  grain  with  an 
outside  diameter  equal  to  that  of  the  forward.  If  this  diameter 
proved  to  be  unsatisfactory,  larger  or  smaller  diameters  were  used. 

For  forward- chamber  designs  with  high  length-to- diameter  ratios,  the 
aft-chamber  subroutine  considered  the  feasibility  of  placing  the  aft 
chamber  around  the  forward  chamber,  rather  than  in  tandem,  to 
reduce  the  over-all  length-to-diameter  ratio  for  the  dual-chamber 
system.  Figure  4  shows  some  typical  configurations  for  chambers 
in' tandem  and  with  one  chamber'  inside  the  other. 

The  subroutine  assumed  that  the  surface-time  history  of  the 
selected  grain  configuration  was  neutral.  This  assumption  is  true 
for  end-burning  grains  and  can  be  made  approximately  so  for  the 
star  and  wagon-wheel  configurations.  Although  most  of  the  cylindrically 
perforated  designs  would  not  be  neutral,  this  can  be  corrected  in  a 
detailed  design  by  using  slotted-tube  configurations  of  approximately 
the  same  diameter  and  length. 

4.  WEIGHT  SUBROUTINES 
a.  General 

As  mentioned  in  Section  I,  the  weight  calculation  sub¬ 
routine  was  divided  into  two  parts  (A  and  B)  to  facilitate  computa¬ 
tion  in  NorthrofrCarolina's  IBM  1620  computer.  The  subroutine 
calculated  the  weight  and  dimensions  of  the  forward  and  aft 
motor  case,  end  closures,  insulation,  nozzle,  igniters,  and 
the  control  valve.  The  weights  and  dimensions  were  selectively 
calculated,  depending  on  geometrical  constraints  such  as  the 
forward-  and  aft-chamber  grain  configurations  and  the 
materials  used  for  each  individual  component.  The  generalized 
flow  diagrams  for  the  two  weight  subroutines  are  shown  in 
Figures  5  and  6,  respectively.  The  subroutine  program  listings 
are  given  in  Appendix  A. 
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Figure  6  -  Weight  Subroutine  B  Logic  Diagram 


-20- 


AFRPL-TR-65-209,  Vol  II 


b.  Chamber  Weights 


Chamber  weights  were  calculated  by  multiplying  the 
chamber  surface  area  by  the  thickness  of  the  particular  section 
(based  on  the  operating  pressure  times  material  density.  The 
weight  of  the  head  closures  was  calculated  from 


W 


ellipsoid  =  'ell'”'  CrS  +  &><ln  ]  - 


1  +  e , 


where 


PRSF 

*.u  '  753  •  and 


(7) 


e  = 


The  cylindical  chamber  weight  was  calculated  from 


Wcyl  ■  *c,l  Pa,R’ 


(8) 


where 


_PRSF 
cyl  CT 


The  symbols  used  in  Equations  7  and  8  are  defined  below: 


W 


W 


ellipsoid 

cyl 


ell 

fccyl 


P 

R 

b 


e 

a 

SF 

P 


total  head  closure  weight  (lb^) 
total  cylindrical  chamber  weight  (lb^) 
thickness  of  1*.  ad  closure  (in.  ) 
thickness  of  cylindrical  chamber  (in.  ) 
material  density  (lb^/in.  3) 
major  case  radius  (in.  ) 
minor  ellipse  radius  (in.  ) 
ellipse  eccentricity  ratio 
material  yield  strength  (lb^/in.  *) 
safety  factor 

operating  pressure  (psia) 
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£.  Valve  and  Power  Supply  Weights 

Valve  and  power  supply  weights  were  determined  by 
a  parametric  derivation  of  each  component's  contribution  to 
total  control  weight  as  a  function  of  valve  size.  The  weight 
equations  used  for  each  component  are  listed  below: 

Component  Weight  Equation 


1 


Intermediate  plate 

W  =  0.  00939  Dt3  (AP)8 

Intermediate  plate 

insulation 

W  =  0.  00  25  D  a  t, 

t  D 

(t,  >  9  sec) 

b 

w  =  0.  0225  Dt8 

(t,  <  9  sec) 
b 

Housing  insulation 

W  =  0.  00  233  D  3 1, 
t  D 

(tfe  —  9  sec) 

W  =  0  .  0  2097  D  3 

V 

(t,  <  9  sec) 
b 

Tube  insulation 

W  =  0.  0053  D  3  t, 
t  b 

(tfe>  9  sec) 

W  =  0.  0476  Dt3 

(tfe  <  9  sec) 

Piston  insulation 

W  =  0.  00162  D  3  t. 

t  b 

(t,_>  9  sec) 
b 

W  =  0.  0146  Dta 

(tfe  <  9  sec) 

Pintle 

W  =  0.  0282 

(Dfc  >  1.  07 

and 

L  =  f(Dt) 

W  =  0.  0302  Dt3 

(Dt<  1.07) 

Power  pack 

W  =  0.  33  D  3 

l 

Servovalves 

w  =  0.  59  Dt3  +  0.  9 

Piston 

W  =  0.  0  22  D  *  +  0. 

00676  Dt3 

(Dt>  2.67) 

W  =  0.  0403  Dt3 

(D  <  2.67) 

l 
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Component 
K^d  of  piston 


Cylinder 


Duct 


Struts 

Housing 


Weight  Equation 

W  *  0.  022  0;  0000581  d/  /P 


(Dti  0.68) 

W  =  0.  22  D  *  +  0.  0000294  D  *  /P 

(D  <  0.  68) 

X 

W  =  0.  0297  D  *  (D  —  2.9) 

»  V 

W  =  0.  00792  Dt*  +  0.  027  Dt* 

(Dt  <  2.  9) 

W  =  3D  a  =  0. 00563  Dt* 

(DA  *  10.7) 
t 


W  =  CDfca  *  0.060  2  Dta 

(D.  <  10.  7) 

X 

a 

W  *  0.  000301  DtaP* 

W  =  0.  0966  Dt*  (Dt  >  2.  38) 
W  =  0.  0934  Dta  +  0.  575  Dta 

(D  <  2.  38) 


where 

D  =  /4  V .  (input), 

*  FcT 

t,  =  burn  time, 

D 

AP  =  forward- chamber  pressure  -  aft-chamber 
pressure  (maximum), 

C  =  constant  input,  and 

yf  =  final  free  volume  of  forward  chamber. 
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d.  Insulation  Weights 

Insulation  weight  was  determined  as  a  function  of  exposed 
insulation  area,  flame  temperature,  burning  time,  and  insulation 
density,  by  means  of  the  following  equation: 

W.  =  A..  p  t.  , 
ins  ins  ins 


where 


T  +  460) 

t.  =  0.  007  t,  I ,  ~--v- 
ins  b  \  5806 


+  0.  05, 


A.  =  exposed  insulation  surface  area  (in.  ) 
ins 

p  =  insulation  density  (lb  /in. 3) 

m 

t,  =?  maximum  duration  (sec),  and 

D 


Tc  =  chamber  flame  temperature  (  F). 


e.  Nozzle  Weights 

Nozzle  weight  was  determined  by  the  following  empirically 
derived  equation,  which  is  dependent  upon  the  yield  strength  of 
the  material  used  in  each  component  of  the  nozzle,  flame  tem¬ 
perature,  expansion  ratio,  throat  area,  chamber  pressure, 
maximum  burn  time,  and  propellant  C*: 


W 

noz 


■  |104-^fe5<MAt)1,,pe  +l.53(At, 

max 


0*9 


°.  00268  tb(1fr^1"’(Pc 


)°*  8  +  0.  1  -I 


max 


10® 


(0.0865)  (€-4)  A., 

I 

where 

€  =  expansion  ratio, 

A^  =  throat  area  (in.8 ) 
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Pc  =  maximum  chamber  pressure  (lb  /in.  3  ) 

max 

=  maximum  burn  time  (sec),  and 

* 

C  =  propellant  characteristic  velocity  (fps). 

£.  Subroutine  Output  Data 


Output  data  from  these  subroutines  was  very  complete. 
An  attempt  was  made  to  summarize  all  motor  characteristics, 
dimensions,  and  weights  in  one  tabulation,  shown  in  Table  III. 
Output  data  are  coded  by  a  line  heading  that  increases  from 
100  to  1800.  ~ 
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1700  Nozzle  Weight  Nozzle  Structure  Nozzle  Throat  Nozzle  Insulation  Expansion 

Weight  Weight  Weight  R»tio 

[800  Total  Motor  Maximum  Motor  Maximum  Motor  Maas  Fraction  Ideal  Velocity 

Weight  Diameter  Length 
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SECTION  m  -  DISCUSSION  OF  RESULTS 


1.  GENERAL 

The  over-all  results  of  the  parametric  study  are  presented 
and  discussed  in  this  section.  The  data  from  the  computer  runs 
have  been  plotted  in  semi-logarithmic  form  to  facilitate  analysis. 
The  data  from  these  basic  plots,  which  are  presented  sequentially 
in  Appendix  B  by  phase  and  variable,  have  been  cross-plotted  and 
presented  in  this  section  to  highlight  the  trends  observed. 

This  section  is  organized  by  study  phase,  with  the  effect  of 
each  independent  variable  discussed  separately. 

2.  PHASE  I 

a.  General 

In  Phase  I,  the  effect  of  specific  impulse,  thrust  ratio, 
number  of  starts,  total  impulse,  and  minimum  thrust  on 
mass  fraction  and  delta  velocity  was  studied.  The  effect  of 
each  of  these  variables  is  discussed  below. 

b.  Specific  Impulse 

The  effects  of  specific  impulse  on  mass  fraction  and 
delta  velocity,  for  a  thrust  ratio  of  1  and  10  starts,  are 
shown  in  Figures  7  and  8,  respectively.  As  shown  in  Figure 
7,  mass  fraction  decreases  slightly  with  increasing  specific 
impulse.  However,  delta  velocity  increases  as  specific 
impulse  is  increased.  Curves  for  higher  thrust  ratios  and 
other  than  10  starts  are  similar,  but  the  mass  fraction  and 
delta  velocity  levels  differ.  (Confidential) 

c.  Thrust  Ratio 

In  Figures  9  and  10,  the  effects  of  thrust  ratio  (that 
is,  maximum  thrust/minimum  thrust)  on  mass  fraction  and 
delta  velocity  are  shown,  for  10  starts  and  a  specific  impulse 
of  280  lb.-sec/lb  .  For  thrust  ratios  up  to  5,  the  penalty  on 

T  tYl 
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Figure  7  -  Effect  of  Specific  Impulse  on  Mass  Fraction 
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Figure  9  -  Effect  of  Thrust  Ratio  on  Mass  Fraction 
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mass  fraction  and  delta  velocity  is  small.  This  penalty  becomes 
more  pronounced,  however,  for  a  thrust  ratio  of  20,  particularly 
for  the  larger  motors.  For  example,  the  mass  fraction  is 
decreased  2.  7  percent  at  a  20-to-l  thrust  ratio  for  a  10,  000- 
lb-sec  motor,  compared  to  a  decrease  of  4.  5  percent  for  a 
1,  000,  000-lb-sec  motor.  This  penalty  is  due  to  the  effect  of  the 
higher  chamber  pressure,  required  for  the  20-to-l  thrust  ratio, *. 
on  case  thickness,  resulting  in  a  greater  weight  penalty  in  the 
large-diameter  (high  total  impulse)  motors.  However,  the 
20-to-l  throttling  range  imposes  very  little  penalty  over  that 
of  a  stop- restart  motor  alone  (thrust  ratio  of  1)  at  the  lower 
total  impulse  levels.  (Confidential) 

d.  Number  of  Starts 

The  effects  of  the  number  of  starts,  for  a  restartable  motor, 
on  mass  fraction  and  delta  velocity  are  shown  in  Figures  11 
and  12,  respectively.  The  differences  represent  the  weight 
penalty  imposed  by  the  multiple  pyrogen  units  for  restart 
capability.  The  weight  of  a  valve  for  termination  was  included 
for  all  motor  designs,  including  the  motor  with  one  start.  The 
mass  fraction  penalty  for  20  Starts,  compared  to  one  start,  is 
1.  8  percent  for  the  10,  000-lb-sec  motors,  and  1.  2  percent  for 
the  1,  000,  000-lb-sec  motors.  (Confidential) 

e.  Total  Impulse  and  Minimum  Thrust 


The  effects  of  total  impulse  on  mass  fraction  and  delta 
velocity  are  shown  in  Figures  13  and  14,  respectively,  for 
motors  with  thrust  ratios  of  1,  5,  and  20  to  1.  The  mass 
fraction  for  each  total  impulse  level  corresponds  to  that  at  the 
optimum  minimum  thrust  level,  derived  from  Figure  B-5.  Figure  15 
is  a  contour  map  of  mass  fraction  as  a  function  of  total  impulse 
and  minimum  thrust  for  a  thrust  ratio  of  one.  A  similar  plot 
for  a  thrust  ratio  of  20  is  given  in  Figure  16.  The  dotted 
constant  mass  fraction  lines  in  Figures  15  and  16  are  extrapo¬ 
lated  values.  .  These  plots  give  the  optimum  minimum  thrust 
levels  corresponding  to  various  total  impulse  level  motors.  (Confidential) 

Figure  17  shows  the  effects  of  minimum  thrust  and  total 
impulse  on  motor  diameter.  The  dotted  lines  connect  the  values 
for  constant  total  impulse  designs,  and  the  solid  lines  connect 
values  for  motor  designs  with  similar  grain  configurations. 
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Figure  11  -  Effect  of  Number  of  Starts  on  Mass  Fraction 
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Figure  12  -  Effect  of  Number  of  Starts  on  Delta  Velocity 
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Figure  13  -  Mass  Fraction  as  a  Function  of  Total  Impulse 
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Figure  14  -  Delta  Velocity  as  a  Function  of  Total  Impulse 
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Figure  15  -  Mass  Fraction  Contour  Map  for  a  Thrust  Ratio  of  1 
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Figure  17  -  Motor  Diameter  as  a  Function  of  Minimum'  Thrust 
and  Total  Impulse 
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A  minimum  diameter  is  obtained  for  each  total  impulse  level 
motor  at  different  minimum  thrust  levels.  Note  that  the  optimum 
minimum  thrust  level  increases  with  increasing  total  impulse. 

In  Figure  18,  the  effects  of  total  impulse  and  minimum 
thrust  on  motor  length  are  shown;  the  dotted  lines  connect  designs 
with  similar  grain  configurations.  This  graph  shows  that  total 
impulse  affects  grain  length  much  less  than  minimum  thrust. 

This  is  also  borne  out  in  Figure  19  in  which  the  band  of  motor 
length  is  plotted  as  a  function  of  minimum  thrust. 

3.  PHASE  II 

a.  General  — 

The  results  of  Phase  II,  in  which  case  material,  burning 
rate  constants,  pressure  exponents,  and  weight  flow  ratios 
were  investigated,  are  presented  in  the  following  paragraphs. 

b.  Case  Material 

Three  case  materials,  fiberglass,  titanium,  and  steel, 
were  investigated  in  Phase  Ila.  The  effect  of  the  mass 
fraction  ratio  of  fiberglass  to  steel  motors  and  titanium  to 
steel  motors  as  a  function  of  minimum  thrust  is  shown  in 
Figure  20  for  100,  000-lb-sec  total  impulse  motors.  The 
effect  of  delta  velocity  ratio  for  these  same  conditions  is 
shown  in  Figure  21.  The  advantages  of  fiberglass  and  ' 
titanium  over  steel  become  more  significant  at  the  higher 
thrust  levels.  (Confidential) 

The  effects  of  case  material  on  mass  fraction  and  delta 
velocity  as  a  function  of  thrust  ratio  are  presented  in  Figures 
22  and  23,  respectively.  As  thrust  ratio  increases,  which 
corresponds  to  higher  chamber  pressures,  the  effect  of  the 
three  materials  on  mass  fraction  and  delta  velocity  becomes 
more  pronounced.  (Confidential) 

c.  Burning-Rate  Constants 

The  effects  of  forward-  and  aft-propellant  burning-rate 
constants  on  optimum  mass  fraction  and  optimum  delta  velocity 
are  shown  in  Figures  24  and  25,  respectively,  for  100,  000-lb-sec 
motors.  For  this  motor  size  the  highest  forward  rate  constant 
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Figure  18  -  Motor  Length  as  a  Function  of  Minimum  Thrust 
and  Total  Impulse 
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Figure  20  -  Effect  of  Case  Material  on  Mass  Fraction 
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Figure  23  -  Delta  Velocity  as  a  Function  of  Thrust  Ratio 
for  Three  Case  Materials 
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and  lowest  aft  rate  constant  investigated  produced  the  highest 
mass  fraction  and  boost  velocity.  These  constants  gave  motor 
designs  with  smaller  diameters  than  others,  as  shown  in 
Appendix  B.  The  higher  forward  rate  constant  permitted  an 
end-burning  forward-grain  configuration  to  be  used,  while  the 
lower  aft  rate  constant  allowed  a  star  aft-grain  configuration 
to  be  employed.  These  optimum  mass  fraction  and  delta  velocity 
values  were  calculated  at  a  minimum  thrust  of  500  lb^.  At.  lower 
minimum  thrust  values,  lower  forward  rate  constants  were  more 
advantageous.  (Confidential) 


d.  Pressure  Exponents 


As  expected,  the  higher  burning- rate  pressure  exponents 
give  higher  mass  fractions  and  boost  velocities  for  throttleable 
motors,  as  shown  in  Figures  26  and  2P .  At  a  thrust  ratio  of  20. 
the  effects  of  forward  and  aft  exponents  are  very  pronounced. 
For  example,  as  the  forward  exponent  increases  from  0.  60  to 
0.  90  and  the  aft  exponent  changes  from  0.  80  to  1.  00  (Figure  26), 
mass  fraction  increases  from  0.  52  to  0.  81.  For  a  thrust  ratio 


of  5  to  1,  on  the  other  hand,  the  effect  of  the  exponents  is  very 
slight.  (Confidential)——- - 


Aft-to -Forward  Weight  Flow  Ratio 


The  effects  of  ait  to  forward  weight  flow  ratio  (0)  on  mass 
fraction  and  delta  velocity  are  shown  in  Figures  28  and  29, 
respectively.  The  higher  0  values  are  advantageous  because, 
at  these  conditions,  the  aft  chamber  operates  at  a  lower 
pressure  than  the  forward  chamber  and  has  a  higher  aft  propellant 
density.  A  0  of  4  produces  less  than  a  4-percent  increase  in 
mass  fraction  over  a  0  of  2,  but  provides  a  7.  5-percent  increase 
in  delta  velocity  (Confidential) 
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Figure  26  -  Effect  of  Burning-Rate  Pressure  Exponents 
on  Mass  Fraction 
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Figure  27  -  Effect  of  Burning-Rate  Pressure  Exponents 
on  Delta  Velocity 
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Figure  29  -  Effect  of  0  on  Delta  Velocity 
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SECTION  IV  -  SUMMARY 


Although  the  performance  values  for  mass  fraction,  delta  velocity, 
and  envelope  generated  in  this  study  are  not  intended  to  represent  those 
that  would  be  obtained  from  detailed  designs,  the  values  agree  fairly  well 
with  those  previously  calculated  in  detail  for  individual  motors. 

It  was  found  that,  by  increasing  specific  impulse,  mass  fraction  is 
reduced  but  delta  velocity  and  total  motor  weight  are  increased.  For 
10,  000-lb-sec  motors,  an  increase  in  specific  impulse  from  265  to  300 
lb-sec/lb  increased  delta  velocity  by  7  percent,  whereas  for  1,  000,  000- 
lb-sec  motors,  the  increase  was  12  percent.  (Confidential) 

The  mass  fraction  penalty  imposed  by  a  20-to-l  thrust  modulation 
was  2.  5  percent  for  low-impulse  (10,  000  lb-sec),  low-thrust  (100  lb^) 
motors.  This  penalty  increased  to  above  4  percent  as  impulse  increased 
to  1,  000,  000  lb-sec  and  thrust  to  4,  000  lb^.  These  values  were  based  on 
the  characteristics  of  available  propellants.  Improved  propellant  pressure 
exponents  would  reduce  the  penalties.  (Confidential) 

For  motors  having  stop-restart  capability  (that  is,  an  on-off  valve), 
the  number  of  restarts  available  had  little  effect  on  mass  fraction  and 
delta  velocity.  The  mass  fraction  penalty  for  additional  igniters  was  less 
than  0.  1  percent  per  restart. 

At  each  impulse  level,  there  was -an  optimum  thrust  range  for 
throttleable  motors  which  was  close  to  the  optimum  thrust  level  for  non¬ 
throttling  motors.  If  the  thrust  was  increased  or  decreased  beyond  this 
range,  both  mass  fraction  and  delta  velocity  decreased. 

Fiberglass  and,  to  a  lesser  degree,  titanium  were  found  tc  be  very 
advantageous  for  motor  case  materials,  compared  to  steel,  for  motors 
with  high  thrust  levels  and  high  throttling  requirements. 

The  highest  forward  propellant  burning-rate  constant  and  lowest 
aft  propellant  burning-rate  constant  investigated  gaive-'fhe  highest  mass 
fraction  and  delta  velocity  values  for  the  100,  000-lb-sec  motor  size. 

Mass  fraction  was  more  sensitive  to  forward  and  aft  propellant 
pressure  exponents  than  to  any  other  propellant  characteristics  at  high 
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throttling  ratios.  Aft  propellant  pressure  exponent  significantly  affected 
mass  fraction  for  motors  with  low  forward  exponents.  However,  with  high 
forward  exponents,  the  effect  of  the  aft  exponent  was  much  less.  Conversely, 
forward  exponent  change  is  less  significant  with  high  aft  exponents.  (Confidential) 

Increasing  the  aft-to-forward  grain  weight  flow  ratio  from  2  to  4,  for 
the  100,  000-lb-sec  motor  size,  increased  mass  fraction  from  0.  80  to  0.  83. 

(Confidential) 
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APPENDIX  A  -  COMPUTER  PROGRAMS 


A  printout  of  the  entire  computer  program  is  presented  on  the  following 
pages.  The  program  is  arranged  by  subroutine  in  the  sequence  shown  in 
Figure  1.  The  individual  subroutines  are  shown  in  more  detail  in  Figures 
2,  3,  5,  an,d  6,  respectively...  --  — 
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-6600  C 


-6600  c 


-6600 

-6600 

-6600 

-6600 

-6600 

-6600 

-6600 


-  ■> 


*600 


-6600  C 


-6600  C 
-6600  C 


-6600  C 
-6600  C 


-6600  C 


-6600 


-6600 


-6600 
/ 


:>600 


-6600 

-6600 

-6600 

-6600 

-6600 


-6600  C 


-5600 

-6600 

-6600 

-6600 

-6600 

-6600 

-6600 

-6600 

-6600 

-6600 


-6600  C 
-6500  C 


-6600 

-6600 

-665? 

-6704 

-6754 

-6708 

-6954 

-6906 

-6950 

-7006 

-7060 


-7199 

-7300 

-73?? 

-7378 

-7436 

-7540 

-756? 


C5R  STEADY  STATE  BALLISTIC  PROGRAM,  FOR  PARAMETRIC  STUDY 
DR.  ALLEY  M.K .  OPPRECHT  AUGUST  1^64 
REVISION  BY  V/.  OSTERHOUT  OCTOBER  1,1964 


6600  C 
C 


S  CODE  1  COOL  FOR  OPERATION  OF  STEADY  STATE  PROGRAM 

5  ?  ”  C0^E  FOR  OPERATION  of  GRAIN  CONFIGURATION  PROGRAM 

S  RHOF  »  FV/D.  PROP.  DENSITY 

S  ALPF  =  FV/T).  PROP.  BURN  RATE  CONSTANT . . - _ 

S  ENPF  «  FV/D ,-  PR-OP-;-f 7(P0tttrrT - 

S  RHOA  *  AFT.  PRO?.  DENSITY 
S  ALPA  »  AFT  PRO?,  BURN  RATE  CONSTANT 
S  EMPA  =  AFT  PRO?.  EXPONENT 
S  CDF  «  FV/D.  PROP,  Cl) 

S  DELPF  =  DELTA  FV/D  PRESSURE 


FMN  =»  MINIMUM  THRUST 
FMX  «  MAXIMUM  THRUST 
PFMU  r,  FV/D  CHAMBER  MINIMUM  PRESSURE 
THEPR  «  THETA  PRIME 
PAMN  -AFT  CHAMBER  MINIMUM  PRESSURE 
XI PR  *  TOTAL  IMPULSE  AT  THETA  PRIME 
S  WF  =  FV/D  PROP  V/EIGHT 
S  ASF  »  FV/D  PROP  AVERAGE  SURFACE  AREA 
S  PFMX  ■  FV/I)  CHAMRER  MAXIMUM  PRESSURE 
S  ASA  -AFT  PROP  AVERAGE  SURFACE  AREA 


IS 


BLANK 
BLA I  !K 
BLAI/K- 
BLANK 
BLANK 


?6)  IS  BLANK 

27)  IS  ATA  =*  AFT  THROAT  AREA 


?8 
29 
-3°, 

C(?1 ) 
C(3?) 
C(33) 

Cf  34) 
C(35) 

19  FORMAT 
19  FORMAT 
70  FORMAT 
73  FORMAT 
24  FORMAT 
?5  FORMAT 
26  FORMAT 


S  V/EBF 
S  WE  BA 
S  WA  « 
S  TBMX 
S  PA  MX 
S  CD  « 
S  ATF 


77  FORMAT 
79  FORMAT 


51  FORMAT 
103  FORMAT 


FV/I)  WEB  THICKNESS 
-AFT  WER  THICKNESS 
AFT  PROP  V/EIGHT 

-  MAXIMUM  V/EB  TIME 

-  MAXIMUM  AFT  PRESSURE 
COO 

>  MAXIMUM  OPERATING  FV/D  THROAT  AREA 
S  BLANK 

1 1HTHETA  PRIME, F14,?) 

11HAFT  SURFACE, FI 4. ?) 

10HAFT  THROAT, F16. 3) 

7HAFT  V/EB, F19. 7) 

13HAFT  PROP.  V/T.  ,F17.7) 

1  1HFV/D  SURFACE, FI 4. 7) 

7HFWD  V/EB, FI .9. 5) 

13HFV/D  PROP.  V/T., FI 7. 7) 

17HMAX  V/EB  TIME,F13.7l 
7HFWD  V/EB, FI 9. 7.) 

(30H 

P  AFT 


30  F0RMAT(/36H  P  FV/D 

31  FORMAT  (/) 

35  FORMAT  ( / 1 3HPA  GREATER  PF) 


THETA 


,F70.6,»5) 

THRUST) 


43  FORMAT!  F9.1.F9.1  F9.?,F9.1.F9.7,F9.0,F9.?,F9.4V 
57  F0RMAT(36X,?6H  WT  AFT  IMPULSE  WEB  it  FT  AT  FWD 

A  4  •  I  % 


101  FORMAT(  14) 

-?0. 


10?  FORMAT!  3F70.10) 
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-7S9A 
-innn 
■  7VV- 

-7r»'>U 
-7S-U 
-7S3L 
-7  SOS 
-7',n 


C  SWITCH  1  OFF  TO  3*»rtCH  C(n  THRU  c(?5) 

C  WITCH  ?  OH  TO  ACCEPT  INPUT  CHANGES  ANYWHERE  I!!  CfO  THRU  C(?5' 

C  SWITCH  3  '!N  STOPS  /LI.  OUTPUT  PRINTING  EXCtPT  SD.W-IARY  LI  fit 

C 

o  i  me  us  i  on  the  f  30 ) ,  X I SP  ( 30 )  ,C •>(  30 )  , C ( 35 ) 
so;  PR  I  ITT  ?1 
READ  10? 

PR  HIT  103 


-ts-to- 

-1V\?  KK=1 

-73S’-  J  FfSEUSE  SWITCH  ?>S17, 

-7.374  sio  pr i riT  31 

-7333  READ  101, N 

-7710  DO  ADD  K»l/I 

-71??  ADD  READ  lO?,THE(K^,XISPfKT,C»i(K) 

-7  (71  DO  SOI  jai.35 

-7T70  READ  103  C{.l) 

-1777  |F(C(.n)SS0,S01,SS0 

-3013  SSO  PR  I  NT  103,  CU) 

-•toss  soi  coiiTimiE 

-nor  so  to  sro 

-3110  Sir  RE/.D  103.X, HP 

-3143  lr(x)51 1 ,5?0,si 1 

-3ror  sii  pr 1 nt  i03,x,  up 

—3?  33  C(NP)~X 

-3?7A  GO  TO  sir 

-3? 9?  sro  PRINT  31 

-3?9A  sn=C{ 13) 

-3303  SS~C(1M 

-3313  S6=C( 15) 

-3330  NF=*0 

-334?  LL=C( 1 ) 

-3373  F1=»C(  11) 

-3390  P1=C( 13)  —  % 

-940?  go  to(io, :i,?i),ll 

-3 A3?  C  10-OPTlOM  I A 

-3A3?  C  1 1=»0PTI  ON  10 

C  10  C(?0)“CO'l»)*C(3)*C(  10)*C(,»)*Ctl3)**C(5)/(C(6)*C17)*C(15)**C(-")) 

Zqtlh  ti  c(Wrf(S)*c(»o)*c(7)*c(i5)«c(»)/(e(s)*c(in)*c(fc)*c(ijj«c(5)) 


-3303  1?  M=*1 

-3910  13  THE  I  **0(14) 


-39?? 

-9934 

-90?  5 

-90S? 

-9070 

-9??5 

-93S3 

-9*v90 

-9S70 

-976? 

-Q73S 

-9910 

-9934 


r^ss  1 

14  IF(THE(K+1)-THEI}1S,16,1& 

15  K»K+1 

CO  TO  1 4 

IS  RAT  1 0= ( THE  I -THE  00  ) / ( THE ( K+1 )-THE (  K  )  ) 

XISPO-X|SPfK)+RATIO*(XISP(K+n-XISPJK)> 

C(33)-CD(K)+RATI  0*f  C«i(K+1  )-CD(K)) 

17  C(7^)i!i'.rtll^5*fC(3)*Cfl«'*C'M*Cfl3)«C(?)V(C(??l*C(tO) 

PRINT  1 9,C( 1 4) 

PRINT  19,C(?0) 

SHW?iWRu>*c(m> 
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-9337 

-990  6 

-991  4 
-9976 
-9934 
-999 4 
J0030 
J0066 
J0176 
J0 1.36 
J0734 
J0782 
J0650 
J0474 
JO  49  9 
J0677 
J0546 
J0570 
J0594 
J0619 
J06'i7 
^694 
J069*>^ 
J0679 
J0774 
JO  3  70 
J0954 
J  107.7 
J 1034 
J1047 
J1134 
J1776 
J1774 
J173? 

J 1 3 1 3 

J  1410 
J1502 
J1550 
J1562 
J1S70 
J1630 
J1666 
J1714 
JI858 
J1918 
J1994 
J7074 
J7074 
J2074 
J214? 
J2210 
J2730 
J2338 
J2394 
J7406 
J2418 
J7430 


PRINT  51,C(?8) 

GO  TO  79 
M=? 

GO  TO  13 

C(  1 7 )=»C ( l6)/(XI  SPO*(  1  ,+C(  14)}) 

C(30)=C(17)*C(14) 

C ( 3 1 )»C( 16)/C( 1 1 ) 

C(?.8)-C(31  )*(C(4)*C(  13)**C(5>) 

C(29)«C(31 j*(C(7)*C( 15)**C(8>) 

C(  18)=*C(  17;/(C(3  j*C(78)  ) 

C(?.0)-C(30)/(C(6)*C(29)) 

PR?NT”20‘c(27)^*C{3>*C(l8)*C(lf)*C(13)**C(5)/(Ct?3)*C(t5)) 
PRINT  1 9 )  C  (  7.0 ) 

PRINT  73,C(?9) 

PRINT  74,C(30) 

PRINT  75,C(18) 

PRINT  26,C(78) 

PRINT  77,C(17) 

PRINT  78,C(31) 

PRINT  30 
PRINT  57 
PRINT  31 

COm=C(?}*C(l8)*C(4)*C(l3)**C(5^ 

'^«13»C  ( 6)*C  (  20  ^*C{  7)*C(1 5  )  **C  ( 8) 
PACAb<CON1+CON3)/(C(77)*C(33)) 

I F(PACAl>CXl3))50,34  34 
PRINT  35 
GO  TO  46 

IFfPACAL-.998*cn5m^36.48 
I  F(  PACAL-1 .007*C  (16)  )3£S4.  49 
C(15)=7.*PACAL-C(15) 

GO  TO  3? 

THECA=*C  0N3  /C  ON  1 
I  F(THECA-.998*C(1 4)  >38,39,37 
IF(THECA-1.00?*C(14))39,39,38  \ 

C(  >0.5*(C(14)+THECA) 

M-3  \ 

GO  1  13  \v 

F-Xl  >i,0*(C0W1*(  1  .+THECA)) 

C(30)«C(17)*THECA  \ 

G(34 )«C0N1/(C(9j*C( 13))  ^ 

C(7.9)-C(28)*C(7)*PACAL**C(8)/(C(4)*C(13)**C(5)) 

Xf»C( 1 7)*( 1 .+THECA )*XI SPO 
GO  T0(509,508),KK 
GO  10(40,40,41 ),LL 

40- 0PT I  OH  I 

41 - OPT I  OH  II 
IF(C(13)-C(19))47.,42,46 
IF(F-C(17))47,42,46 

1F( SENSE  SWITCH  3)513,514 

PRINT  43.C(13), PACAL, THECA, F,C(30), XI ,C(29),C(34) 
IF(NF)1001, 1002, 1001 
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J?A4? 

J  2454 
J?4ftft 
J?50? 
J2510 
J?52? 
J?66ft 
J?9?2 
J7954 
J3159 
J3166 
J3173 
J3786 
J3798 
J3310 
J337? 
J3334 
J5346 
J3353 
J3370 
J3390 
J3M)? 

J  349ft 
J3494 
PROCESS 
START 


1001  F1-F 

P1-PACAL 

C(13)»C(13)+C(10) 
no  TO  44 

4ft  KK=?  — 

C  19 )»  c  1?)/XISP0)-C(6)*C(70)*C(7)*S^ 
C(19)=(C(19)/(C(3)*C(1R)*C(4)))**(1./C(5)) 
^EC4“(C(o)*c(7)*c(7n)*c(3?)**c(n))/(C(3)<fCcit)*c(tp)^c(19)**c(5)) 

508  PRIUT  31 

C(7.qj„s?,C(19>,C(3?)*THECA»C(1?)*C(3°>»X,.Cf?<J).c(34) 

c(3c»r 

C(34). 


(*) 


506 

503 

504 


52 

53 

54 
$5 
S6 


ING 


CC 13) 

Cf  14) 

C(  15) 
PAUSE 
I F(SEHSE 
DO  504  J 
PUNCH 
GO  TO 
END 

COMPLETE 


SWITCH 
1,35 
10?,C(J) 
50?. 


1)502,503 
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-6600  C  PDO  FORTRAN  PROGRAM 
-6600  C 

-6600  C  C{19)  IS  FORWARD  CHAMBER  AVERAGE  SURFACE, SO  IN  ASF 

-6600  C  C(?9)  IS  FORWARD  CHAMBER  WEB  THICKNESS,  III  WECF 

-6600  C  C(?7)  IS  AFT  THROAT  AREA,  SO  |||  ATA 

-6500  C  Cf?0)  IS  AFT  AVERAGE  SURFACE  AREA,  SO  III  ASA 

-6600  C  Cf79)  IS  AFT  WER  THICKNESS,  III  WEBA 

-6600  C 


-6600  1  FORMATION  ,F20.6,I3) 

-6694  2  FORMAT  (F20.6) 

-6716  3  FORMAT  (/22H  FWD  GRAIN  IS  A  CYL  EB) 

-67RO  4  F0RMAT(1?H  FWD  PROP  0D.F10.2) 

-6844  5  F0RMAT(13H  FWD  PROP  LEN,FI0.2) 

-6900  6  FORMAT? 13H  FWD  PROP  L/D,F10.2) 

-6956  7  F0RMAY(?22H  FWD  GRAIN  l§  A  SPH  EB) 

-7030  8  F0RMAT?/27.H  FV/D  GRAIN  IS  A  CYL  CP) 

-7104  9  FORMAT? 18H  FWD  CHAM  VOL  L0AD,F10*4) 

-7170  10  F0RMAT(/?2H  FWD  GRAIN  IS  A  SPH  CP) 

-7744  11  F0RMAT(/24H  FWD  GRAIN  IS  A  CYL  STAR) 

-737?  12  F0RMAT(/24H  FV/D  GRAIN  IS  A  SPH  STAR) 

-7400  13  F0RMAT(/22H  FWD  GRAIN  IS  A  CYL  WW) 

-7674  102  F0RMAT(/19H  AFT  GRAIN  IS  AN  EB) 

-7542  103  F0RMAT(14H  AFT  GRAIN  DIA,F10.3) 

-7600  104  F0RMAT(14H  AFT  GRAIN  LEN.F10.3) 

-7659  105  FORMAT? 19H  AFT  GRAIN  VOL  L0AD.F10.3) 

-7726  106  FORMAT? 15H  AFT  GRAIN  P0RT.F10.3) 

-77^6  107  FORMAT? /IBM  AFT  GRAIN  IS  A  CP) 

-7952  103  F0RMAT?25H  AFT  PORT  TO  THROAT  RATI0.F10.3) 

-7R3?  109  F0RMAT(/?0H  AFT  GRAIN  IS  A  STAR) 

-3002  110  FORMAT? /18H  AFT  GRAIN  IS  A  WW) 

-8068  111  F0RMAT?/39H  AFT  GRAIN  IS  AN  EB  OUTSIDE  OF  FWD  CHAM) 

-8176  112  FORMAT?/38H  AFT  GRAIN  IS  A  CP  OUTSIDE  OF  FV/D  CHAM) 

-323?  500  FORMAT  (30H  ) 

-3366  DIMENSION  C(37) 

-3366  C  FORWARD  CHAMBER  CALCULATIONS  OR  PYROGEN  CALCULATIONS 

-9366  1002  READ  500, tj 

-8390  PUNCH  500, N 

-9414  PRINT  500  N 

-3438  CONTROL  102 

-3450  DIAA-O. 

-3462  XLODA-O. 

-3474  CONTROL  102  __ 

-3436  DO  501  J=»1 ,37 

-8493  501  C(J)»0. 

-3570  DO  21  K«1,35 

-8582  21  READ  2,C(K) 

-8666  J=1 

-8678  L-1 

-8690  D!AF=(C(l8)/.7854)**.5 
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-37*9 

-9790 

-9795 

-3994 

-9966 

-9979 
-390? 
-99?  6 
-9950 

-9959 

-90?  6 

-91 ';6 

-9170 

-919? 

-9,9'v 

-9?  19 

-9?  4? 

-9?54 

-9?  4? 

-9?  99 

-9910 

-93  *‘6 

-9454 

-9550 

-9596 

-9654 

-9730 

-9766 

-9779 

-979^ 

-9799 

-9910 

-99?? 

-9946 

-99  "»0 

-9994 

-9919 

*•99?  6 

-996? 

-9974 

-9999 

J0 106 

J0 14? 

JO?  10 

00?  36 

J0?99 

J0334 

J035S 

JO366 

J0379 

J0390 

J040? 

J0410 

JO  446 

J0554 


XLF«C(?3) 

XLOnF-XLF/DIAF 

IF(XLO!>F-.5)?3,??,?? 

??  PRIHT3 
C(37)=?. 

PR I  f'T\ T)l  AF 
PRH!T5,XI.F 
PRIf!T6,XL0!)F 
GO  TO  45 

?3  IF(XL0DF-.?)?6.?4  ?4 
?4  9|AF=((C(l9)*Cf?91  V.M39 
VLF=.9 
PR  I  NT  7 

C(?7)»3. 

4?.  PR!NT4;niAF 
PRINT9.VLF 
XLODF-1 . 

GO  TO  45 

?5  niAF=C(73)/.?5 
M«1  •— 

?5  DIAFS=D!AF*DIAF 

VLF=*(DIAFS-(IMAF-(?  ,*C(?9)  )  )**?  .  )/DIAFS 
XLF*»C(  l31*C(?.S)/(  #7B54*VLF*DIAFS) 
XLODF=XLF/DIAF 
IF(XLOUF-3.)79,?7,?7 
?7  GO  TO  (?3,3t),M 
?3  I)!AF=*C(?3}/.3 
M=?. 

GO  TO  7.6 
?9  PR  INTO 
C(37)»1. 

J=? 

7?  PRINT4fniAF 
PRIHT5.XLF 
PRINT6,XL0I>F 
PRINT9.VLF 
GO  TO  45 

31  D!AF«C(?3)/.?5 
M»1 

VLFa.35 

3?  XLF»(C(18)*C(?9))/(.7R54*VLF*DIAF*DIAF) 
XL0DF=*XLF/I5!AF 
IF(XL0DF-3. )35,33,33 

33  GO  TO  (34,37), M 

34  M=*? 

DIAF-C(?9)/.2 
VLF-.8 
GO  TO  3? 

35  PR  I  NT 11 
C(37)**1  • 

J-3 

GO  TO  7?. 

37  DIAF-C(2B)/.15 

39  XLF«C(18)*C(28)/(.7054*VLF*DIAF*DIAF) 
XLODF-XLF/DIAF 
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JO  5  90 

JO  65  3 

J0670 

J063? 

J0690 

JO  7  70 

J0906 

J0313 

J034? 

J03S0  C 

J0350 

J0O06 

J093? 

J106? 

J 1 1 9? 
J1706 
117% 
J1373 
J 1 390 
JIM)? 
J14?6 
J1450 
J146? 

J 1  474 
J1530 
J167? 

J1 634 
J1646 
J16S3 
J1666 
J 1 67*3 
J1690 
J1753 
J1830 
J1%6 
J1934 
J1946 
J1932 
J2066 
J70V3 
J7090 
J2114 

J7133 

J2162 

J2186 

J7194 

J2706 

J77.42 

J7310 

J2373 

J246? 

J7554 

J2638 

J27% 

J2782 


IF(XL0DF-3.)40,4l,4l 

40  PR  I  NT 13 
C(37)-1. 

GO  TO  72 

41  GO  TO  (45  42,40)  ,M 

42  DIAF=C(?8)/.1 
M=»3 

VLF-.7 
GO  TO  33 

AFT  GRAIN  DESIGN  IF  C(?)-1. 

49  IF(C(?))51 ,73,51 

73  GO  TO  (81 ,20) ,L 
31  GO  TO  ( ?0 ,502, 50? ), J 

niAF=»(c(i9)*c(?n)?.4i%)**{i./3.) 

VLFs.3 
GO  TO  (20,74,77) .J 

74  IF(C(73)-.3*DI AF)79,75,2.0 

75  PRINT  10 
C( 37 ) — 1 . 

76  PRINT  4,DIAF 
PRINT  9,VLF 
XL0DF=1 . 


507 


IF(C(2) )80,20,90 
77  IF(C(2.8)-.2*DIAF)78, 20,20 
73  PR  I  NT 1?  — 

C(37)-1. 

XL0DF=»1  . 

GO  TO  76 
51  CONTROL  102 
CONTROL  102 

30  IF(XLODF-3.)65,65.63 
65  DIAA-((C(20)+2.*C(27))/.786)**.5 
XLODA=*C  ( 29  )  /D I AA 
IF(XLODA-.5)55,53,53 
53  XLA-C(29) 

APA»?.*C(27) 

VLA-1 ,-APA/(DIAA*D1AA*.7B6) 

PR  I  NT 102 


C(36)-2. 

54  PRIHT103,DIAA 
PR1NT104,XLA 
PR1NT105.VLA 
PR  I NT 106, A PA 
GO  TO  73 

55  DIAA-DIAF 
WFA»C(?9)/DIAA 


56 


WFA-»3)60,56,56 

WFA-,4)90,59,59 


11  \  11 1  n-|  -r  / 

90  APA-(I)IAA-?t*C(?9))**?.*.786 
IF(APA-2.*C(27))59,57.57 
57  VLA»1.-APA/(D1AA*DIAA*,786) 

XLA»C ( 2 9 ) *C  2.0 ) / ( , 786*D IAA*DI AA*VLA ) 
PTA-APA/C(27) 

PRINT107 
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j?bo6 

J2930 
J23S4 
J?379 
J290? 
J2910 
J2946 
J304? 
J3126 
J3234 
J3246 
J3753 
J3270 
J327« 
J330? 
JR  374 
J3Tt66 
J3502 
JR510 

jbgir 

J3636 

J369B 

J'3710 

J371R 

J3730 

J374? 

J3750 

J3942 

J3934 

J3970 

J4102 

J4114 

J4162 

J41 74 

J4186 

J4194 

J4230 

J4393 

J45^2 

J4590 

J4638 

J4650 

J4662 

J4670 

J4670 

J4690 

J4702 

J4714 

J4726 

J4762 

J4774 

J47R6 

J49R2 

J4894 

J4902 


Cf  3SW1 . 

53  PR|MT10?tniAA 
PRIMT104  XLA 
PR IHTIOf.VLA 
PRIUT10R,PTA 
GO  TO  73 
59  APA»2.*C(7.7) 

DIAA«?.*C(29)+(APA/.786)**.5 
VLA»1 .-APA/(D!AA*DIAA*.786) 

XLA**C  ( 29 ) *C ( 20 ) / (  .  706*n  IAA*D1  AA*VLA  ) 

PTA«2.0 
PR  I  NT 107 
C(?6)-l. 

GO  TO  50 
GO  VLA=*.85 

69  APA=(1.-VLA)*D!AA*I)IAA*.736 
IF  (APA-2.*c(27))70,71,71 

70  VLA-VLA-.05 
GO  TO  69 

71  XLA^C ( 20 ) *C ( 29 ) / ( . 736*D I AA*D I A A*VLA ) 

I  F(\7FA-.?)62  ,61 ,61 

61  PRINT109 
0(b6)»1. 

GO  TO  54 

62  PR  INT1 10 
C(3$)«1. 

GO  TO  54 

63  I F(Cf 29)-.6*XLF)64t66I66 
6'1  IF(C(79)-.25£DIAF)6B,65,65 
66  APA-2.*C(27) 

DIAA-.2+((C(20)+APA+.786*DIAF*DIAF)/.736)**.5 

XLA=»C(?9) 

VLA=*C(20)/(C(20)+APA) 

PR  INT1 11 
C(36)»3. 

GO  TO  54 
68  APA»2.*C(27) 

DIAA»2.*C(79)+((APA+0.786*(2.*C(29)+DIAF+.2)**7.)/.786)**.5 

C0NI».786*DIAA*DIAA-.786*(IHAF+O.2)**2.-APA 

VLA-CON I / ( CON l+APA ) 

XLA«C  (  29  )  *C  ( 2.0  )  /CON  I 
PR  INT1 1? 

C(36)«4. 

GO  TO  54 

C  PUNCH  OUTPUT  UNLESS  SENSE  SWITCH  1  IS  ON 
20  I F( SENSE  SWITCH  1)1002,1003 
1003  C(22.)-DIAF 
C(?3)-0IAA 
C(24)«XL0DF 

C(?5)  -  XLA/DIAA  X 

C(?6)«C(?) 

DO  1001  K-1.37 
1001  PUNCH  1 ,C(K) ,K 
PAUSE 

GO  TO  1002 
END 


PROCESSING  COMPLETE 
START 
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-6600  C 
-6600  C 
-6600 
-6600 
-6612 
-6703 
-6744 
-6730 
-6816 
-6884 
-6920 
-6928 
-6964 
-6976 
-6988 
-7000 
-7012 
-7024 
-7036 
-7043 
-70  60 
-707? 

-7084 

-7086 

-7109 

-71?0 

-713? 

-7169 

-7130 

-7192 

-7204 

-7216 

-7223 

-7264 

-7276 

-7233 

-7300 

-7312 

-7324 

-7336 

-7343 

-7360 

-7372 

-7334 

-7396 

-7403 

-7420 

-7437 

-7444 

-7456 

-7463 

-7430 

-7492 

-7504 

-7516 


CSR  COMBINED  WEIGHTS  ANALYSIS  PROGRAM 
PROGRAM  5017A  24  NOV  64  M  A  TODD 

DIMENSION  C(70) 

DO  250  J-1,70 
250  C(J)«0 
10  READ  200.X, I 
PUNCH  200, X, I 
201  READ  200, X, I 

I F(X-9999999. >204,17,1 7 
204  C(I)«X 
GO  TO  201 
17  MCAG=C(36) 

RHOF  »  C(3) 

ALPF  «  C(4) 

ENPF  »  C(5) 

RH0A»C(6) 

ALPA»C(7) 

ENPA»C(8) 

PAMN«C( 15) 

XI PR  -  C( 16) 

WF=C( 17) 

ASF  -  C( 18) 

PFMX«C( 19) 

ASA=C(?0) 

DIAF»C(22) 

niAA=C(23) 

CHBMU  =»  C(?6)+1 . 

ATA=*C(?7) 

WA*C(30) 

TBMX-C(31 ) 

PAMX-C(37) 

ATF=»C(34)  — 

MCFG«C(37) 

R0CHF=»C(33) 

ERATA=»C(39) 

SDCHA»C(40) 

R0CHA**C(41  ) 

R0INF«C(42) 

R0INA»C(43) 

CD~C(44) 

WMIG»C(45) 

GSIG»C(46) 

FLTPI«C(47) 

CPF-C(48) 

THC0F»C(49) 

WMFG-C(50) 

TEMF»C(51) 

C0NIG-C(52> 

XMFIG*»C(53) 

ENIG-C(54> 

SFIG-C(55) 

SDNTH»C(56> 

SDNST-C(57) 

SDNXN-C(53) 

EXPRN-C(59) 
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-•»5?i 

-7540 

-799'i 

-7979 

-7911 

-7900 

-791? 

-7910 

-790? 

-7704 

-771? 

-7714 

-7957 

-7194 

-7979 

-7984 

-7996 

-799?. 

-90?8 

-911? 

-9194 

-8797 

-RA?0 

-8916 

-99?  4 

-99?0 

-1699 

-9796 

-1760 

-9109 

-1969 

-91?  9 

-1996 

-9070 

-9091 

-919? 

-9700 

— 93?0 

-9988 

-949? 

-9900 

-9900 

-996.9 

-965? 

-9774 

-979? 

-9960 

J0056 

J0064 

J0 160 

J0708 

JO?  76 
J0300 
J0349 


TKHAaCf 601 
THC0A*»Cf 6? 1 
ERATF«cMl) 

SDCHF-Cf6  :) 

PRC0!l«Cf  691 
CPA«C(661 
WMAG  =»C(67) 

IF(C(?5)-1 .)31 ,31 ,3? 

31  EL0!)A®0 . 

erata=»i  .0000001 
GO  TO  33 

3?  EL01)A=C(?5)-(ERATA-1 . V?. 

33  IF(C(?4)-i.)34,?4,39 

34  EL0I)F«0. 

,ERATF=1 .0000001 
GO  TO  39 

35  EL0DF=C(?4)-(EHATF-1 .)/?. 

39  RI)CHF=I)I  AF/7 . 

RDCHA-DIAA/7. 

4  RAI)=.(  1  .-fl  ./ERATA  )**?.)**.  5 
V/CX  -  (fRT>CHA/ERATA)**7./7.*RAR) 

1F(ERATA-7.5)  46.46,45 

46  AFCL  -  (WCX*L0G((1.-^AD)/(1.-RAD))4RDCHA**?.)*3.14159 

TFCL-(:-‘AMX*RDCHA  )/( 1  .M4*SDCHA*R0CHA) 

GO  TO  47 

45  TFCL=RPCHA*( (3.*PAMX)/4.*CI>CHA*R0CHA)**.5 
AFCL  «  3.14159*(RT»CHA)**7. 

47  I Ff TFCL- .0? 9 1 41  44 , 44 

43  TFCL=.C?5 

44  WFCL=»A  FCL*TFC!. ’‘ROCHA 

ACVL  «  1?.566*RI)CHA*RT»CHA*E101)A 
TC  Vi  -  ( :>A  MX*RT)CHA  /  SDCHA  »‘F*  CC  &  ) 

|FfTCVL-.0?5)M ,4?,4? 

41  TCVL®.0?5 
4?  WCVL»ACYL*TCYL  •‘ROCHA 

WACL=»f  AFCL-?  ,*VATA  )*TFCL*ROCHA 
WCH  BT«V/FCL-*tfA.C  L+Y/CYL 

T 1  iJS«  ,007*T3HX*  (  (  (TF.HA+460 .  )/5806.  )**'■ .  )+.05 
V/IHFC=AFCL*T  I  I1S*R0I  f!A 
Wl  NAC=(AFCL-?  ,*ATA  )*T1  ?!5*R01  HA 
W I HCY=ACYL*T I NS*RO 1  HA 

IF(CHR?III-?.1?09,8,9 
9  F  A7)®(  1 1  ,/ERATF}**? .  )**.5 
WFCX«( (RDCHF/ERATF)**?.  ./7.*F  AD) 

IF(ERATF-?.5)  460,460,450 

460  FAFCL=(WFCX*LOG{ ( 1  ,+F  AI))/(1.-F  AD))+RPCRF**?.)*3.14159 
FTFCL®(PFMX*RDCHF  )/( 1 ,414*SDCHF*R0CHF) 

GO  TO  470 

450  FTFCL«R1)CHF*( (3 .*PFMX)/4.*SDCHF*ROCHF)**.5 
FAFCL  -  3 . 1 41 59*(RDCHF)**7 . 

470  1 F( FTFCL- ,0?5 ) 430 , 440,440 

430  FTFCL®. 075 

MiO  FWFCL»FAFCL*FTFCL*ROCHF 

FACYL  ®  1?  .566*RDCHF*RI)CHF*EL0I)F 
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JO  40  9 
J049? 
J0560 
J0534 
J063? 
J0716 
J0764 
J0834 
J093? 

J 1 01 6 

J1064 

J1064 

J1144 

J 1 1 44 

J1130 

J1139 

J1138 

J1?36 

J1?44 

J 1 2  44 

J1292 

J1292 

J1376 

J1376 

J1412 

J 1 420 

J1420 

J1468 

J1476 

J1476 

J1620 

J1628 

J1628 

J16R8 

J1 724 

J1760 

J 1 772 

J1734 

J1796 

J1796 

J1380 

J1952 

J7036 

J2120 

J2204 

J2276 

J2360 

J2444 

J2528 

J261? 

J2672 

J2756 

J2816 

J2900 

J2994 


C 

C 


C 


C 

C 

C 


C 


C 


C 


C 


FTCYL«(  PFMX*RI>CHF)  /  ( SDCHF*ROCHF) 
IF(FTCYL-.025T410,420,470 
410  FTCYL-.075 
470  FWCYL«»FACYL*FTCYL*ROCHF 

FWACL»(FAFCL-7.*ATF)*FTFCL*R0CHF 

fchbt-fwfcl+fwacl+fwcvl 

FT  I  NS«  ,007*TBMX*  ( ( ( TEMF+460 . ) /5RO6 .  )**4 . )+,05 
F I NFC=FA FCL*FT INS*R0 I NF 
F I NAC** ( FA  FCL-7 . *ATF ) *FT I NS*R  0 1 NF 
FINCY=FACYL*FT1NS*R0INF 

GO  T0( 210, 2 11, 212) ,MCFG 

FWD  GRAIN  ISA  CYL  OR  SPH  CP,  STAR,  WAGONWHEEL 

210  FWI  NS=*F  I NAC+F I  NFC 
GO  TO  209 

FWD  GRAIN  IS  A  C*L  END  BURNER 

211  FWI  NS=*FINAC+(FI  HCY/7  . ) 

GO  TO  7.09 

FWD  GRAIN  IS  A  SPH  END  BURNER 
7.12  FW I  NS**F  I  NAC+(  F 1  NFC/7  .  ) 


209  GO  T0( 70S, 206, 207, 208) ,MCAG 

AFT  GRAIN  IS  A  CP, STAR, OR  WAGONWHEEL 

205  W I NS-W I NFC+W I NAC 
GO  TO  100 

AFT  GRAIN  IS  AN  END  BURNER 

206  W I  NS«W I  NAC+( W I NCY/7.  •  ) 

GO  TO  100 

AFT  GRAIN  IS  AN  END  BURNER  OUTSIDE  OF  FORWARD  CHAMBER 

207  Wl NS»WI NFC+W I NAC+(WI NCY/2 . )+FI NAC+( FI NCY/2 .) 

GO  TO  100 

AFT  GRAIN  IS  A  CP  OUTSIDE  OF  FORWARD  CHAMBER 
708  W I NS-W I NFC+W I NAC+F I NAC+F I  NFC 
100  WCASE-WCHBT+FCHBT 
WINSL-WINS+FWINS 
WIGT  «  0,0 
WIGP  «  0,0 
WIG  -  0.0 


PUNCH50 1 , SDNST , SDNTH . SDN XN.SDCHA, SDCHF ,  X I  PR 
PUNCH  50? , TEMA ,TEMF, PAMX, PFMX. PAMN 
PUNCH50  3 , ELODA , ELODF , ERATA , ERATF ,D I AA ,D I AF 
PUNCH  504,CHBNU,PRC0N,WA,WF.TBMX,EXPRN 
PUNCH  q05,ATA,ATF,ASA,ASF,R0INA,R0INF 
PUNCH  506,FLTP|,WMAG.WMFG.ALPA.ALPF 
PUNCH  507.SFI G,WMI G.ENIG,gSIG,cONIG,XMFIG 
PUNCH  508,ENPA,ENPF,CPA,CPF-THCOA,THCOF 
PUNCH  509.WI GT.WI GP.WI G.RHOA ,RHOF,CD 
PUNCH  5 1 0 , W I NSL , W I  NFC , W I  MAC , W I NCY  ,W I  NS , MCFG 
PUNCH  511.FI NFC , F I NAC , F I NCY , FWI NS 
PUNCH  512 .WCASE ,WFCL ,WACL .WCYL ,WCHBT ,MCAG 
PUNCH  513,FWFCL.FWACL,FWCYL,FCHBT 
PUNCH  514,TCYL,TINS.RCYL,FTINS,R0CHA,R0CHF 
200  FORMAT (30H  ,F20.6,I3) 

501  FORMAT (4H  100,6F11.0) 
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J3056  50?  F0RMAT(»»H  200,t»F11.U) 

J3119  503  FORMAT ( JtfH  300#6F11.5) 

J3190  504  FORMAT (4M  480,6m  .3) 

J3?4?  505  F0RMAT(4H  500.6F11.3) 

J3304  506  F0RMAT(4H  600,6F11.5) 

j'^366  507  F0RMAT(4H  700, 6F1 1  *4) 

jVtfl  508  FORMAT (4H  8OO/1FII.8) 

j'VvRO  500  FORMAT  (4H  nOO,6F11.5) 

■J355?  510  F0RMAT(4H1000,5F1 1 .3. M 1 ) 

.13614  511  F0RMAT(4H1100,6F11.35 

j’675  51?  F0RMAT(4H1?00,5F11.3  Ml) 

,i*7T8  51*  FORMAT! 4H1 300, 6F1 1  *35 

J3800  514  FORMAT f4H1 400 ,6F1 1  .*•-) 

J3T6?  00  TO  10 

J3H70  EMI) 

PROCESSING  COMPLETE 
START 
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-6600 

-6600 

-6600 

-6600 

-661? 

-6703 

-673? 

-6756 

-6340 

-691? 

-6006 

-7030 

-7164 

-7?36 

-73?0 

-7404 

-7433 

-7577 

-767? 

-771 6 
-7776 
-7360 
-7396 
-7937 
-3000 
-.3000 
-3000 
-3000 
-8000 
-8000 
-31 20 
-3163 
-3704 
-8?.  40 
-8303 
-3344 
-3380 
-3416 
-845? 
-8460 
-8508 
-3556 
-3604 
-9652 
-871? 
-8730 
-391? 
-9q?0 
-905? 
-91?0 
-9156 
-9164 
-9700 
-9?  68 
-935? 


C 

C 


C 

c 

c 

c 

c 

5 


CSR  COMBINED  WEIGHTS  ANALYSIS  PROGRAM 
PROGRAM  50178  ?7  OCT  64  M  A  TODD 

DIMENSION  C ( 3 1 ) 

10  DO  250  J-1,31 
250  C(J)=0 

READ  200.X 
PRINT  200, X 

READ  501 , SRNST , SDNTH , SDNXN , SDCHA , SDCHF , XI  PR 

READ  502 ,TEMA,TEMF,PAMX,PFMX,PAMN 

READ  5O3.EL0DA,ELODF,ERATA,ERATF,DIAA,DIAF 

READ  5  0  A ,  CH  BNIJ ,  PRC  ON ,  WA ,  WF ,  T  BMX ,  E  XPR  N 

READ  505,ATA,ATF,ASA,ASF,ROINA,ROINF 

READ  506,FLTP! ,WMAG,WMFG,ALPA,ALPF 

READ  507 , SF I G ,WM I G , EN I G , GS I G, CON  I G , XMF I P 

READ  503,ENPA,ENPF,CPA,CPF,THC0A,THC0F 

READ  509, Wl  GT.V/I  GP ,WIG,RHOA  .RHOF.CD 

READ  510,WI NSL,W! NFCjWIHAC ,WI NCY  Wl NS,MCFG 

READ  511 ,FINFC,FINAC,FINCYfFWINS 

READ  51? ,  V/C  A  SE  ,WFCL ,  WACL ,  WCYL  ,WCHBT ,  MCA  G 

READ  513,  FV/FCL ,  FWACL ,  FV/CYL ,  FCHBT 

READ  51 4,TCYL .TINS, FTCYL ,FT I  NS, ROCHA ,ROCHF 

RDCHA  =  DIAA/7. 

RDCHF  a  DIAF/2. 

I  F(  CHBf'U— 2 .  )7, 19, 19 
VALVE  WEIGHT  CALCULATION 

VFREE=*C  ( 1 )  ,DTH7«C { ? ) ,  DTH?=*C  ( ? )  tOTH?«C  (  4)  .  PL  I  NS-C(  5  ) 

COVI  N=»C(6)  ,TUBI  N=»C(7/  ,P  •  S  IN*=C(  3)  ,T0INS.C(9),P|  STH«C(  10) 

P I  NT=C(  1 1 )  ,COVE=C(  1 2 ) ,P I  ST=*C(  1 3 )  ,CYL=C(  14)  ,DUCT=C(  15) 
STRUT-C ( 1 6 ) . PL  I  i. '  C ( 1 7 ) . POWPK=C (19), SERVO-C (19), VALVE^C ( 20 ) 
19  C(1)-6.283*(RDCHF«-  ,..)*UL0DF+?./(3.*ERATF)) 


70 


21 


7? 

73 

74 

25 

26 

77 

28 

29 


C(2MATF/.7856)**.5 

C(3)«C(2)**7. 

C(4)«C(2)**3. 

IF(TBMX-9.)20,71 ,2* 

C(5)=.02?5*C(3) 

C ( 5 )«  »02097*C ( 3 ) 

C(7)«.0476*C(3) 

C(3)«.0146*C(3) 

GO  TO  22 

C(5)  »,0025*C(3)*TBMX 
C(6)  =  *002  3  3*C ( 3 ) *T  BMX 
C(7)  =•  .005 3*C  ( 3 )  *T  BMX 
C(3)  =.00162*C(3)*TBMX 
C(9)*»C(5).+C(6)+C(7)+C(8) 

IF(C(7)-. 68)23,24, 24 

C(10M.0?2*C(3W.0000?94*C(3)*(PFMX**.5)) 
GO  TO  75 

C(  10)=(  ,0?.?*C(3)  )+(  ,000058l*C( 4)*(PFMX**,5)) 

IF(C(?).-1  .07)26.77,77 

C(11)a.0302*C(3) 

GO  TO  23 

C(11)=.0?B2*C(4) 

I  F(C(  7  )-?.  .38)29,30,30 

C( 1?)-( .0934*C(3))+( .0575*0(41) 

Gr  TO  31 


A- 14 


AFRPL-TR-65-209,  Vol  II 


-9360 
-9396 
-946': 
-9900 
-9609 
-959? 
-9660 
-9744 
-975? 
-9799 
-9956 
-999? 
-9900 
-9936 
J00?0 
J0116 
J015? 
J0700 
J0370 
J0369 
J0390 
J039? 
JO  40  4 
J0416 
JO  4?  9 
JO  440 
J045? 
J0464 
J0472 
J047? 
J0494 
J0604 
J0616 
J0623 
J0640 
J065? 
J0664 
JO676 
JO688 
JO  700 

JO8O8 
J0928 
J1048 
Jit  68 
J 1 204 
J1240 
J1240 
J1?40 
J  12.40 
J1360 
J1396 
J1437 
J1490 
J1516 
J155? 


C 


C 

C 

C 


30  C(1 2 )=. 0966*0(4) 

31  I F(C( 2 )-2. 67)32,33, 33 

32  C(  13)*».0403*C(3) 

GO  TO  34 

33  C(  13M  .022*0(3) )+( .00676*0(4)) 

34  IF(C(2)-2. 9)35,36. 36 

35  C( 14)«( ,00792*0(3) )+( .027*0(4)) 


GO  TO  37 

36  C(14)«. 0297*0(4) 

37  IF(C(2)-1 .07)39.39,39 

38  C(15)=.060r*C(3) 

GO  TO  40 

39  C( 15)=. 00563*0(4) 

40  C(16)=(.000301*C(4)*(PFMX**.666)) 

C(17)=(. 00939*0  ,'4)*((PFMX-PAMX)**. 5)) 

C( 18)=. 33*0(4) 

C(19)=(.59*CT4))+.9 

C(31)=C(9)+C(10)+C(11)+C(12)+C(13)+C(14)+C(15)+0(16)+C(17) 

o(?o)=»c(3i)+c(in)+c(iP) 


AS^ASF 


RHO=»RHOF 

TEM=>TEMF 

EMP-EMPF 

CP=»CPF 


ALP«ALPF 

THCO-THCOF 

WMG-WHF6 

GO  TO  100 

IGNITER  WEIGHT  CALCULATION 
7  PRINT  70 

C(  1  )p,6 . 2  83*  (RDCHA**3 .  )*(  ELODA+?  ./  ( 3.*ERATA  )  ) 

AS=»ASA 

RH0=RH0A 

TEM-TEMA 

ENP-ENPA 

CP=CPA. 

ALP=»ALPA 

THC0=»THC0A 

WMG-WMAG 

100  C0N1»AS  *RHO  *TEM  *18550./ (0(1 )  *WM  G) 

C0N2-(C0N1*2.*THC0  *ENP  )/(CP  *ALP  *RHO  ) 

PSTAR  «  (C0NIG*C0N2)**(  1 ./( 1  ,+EIIP) ) 

W1GP»(WMIG*PSTAR*C(1)  )/(GSIG*FLTPI*l8550.)*SFIG 

WIG-WIGP/XMFIG 

WIGT»WIG*ENIG 

NOZZLE  WEIGHT  CALCULATION 

WNST-C  ( ?  1 )  .  SDST-C (  2 2  )  .WNSTR-C  ( 23  ) , WNTH-C (  2.4) , SPTH=C  (25) 
WI!THT=»C(?6)  ,WNX.I»C(?7)  ,SDXN*C(7H)  ,WNXlFf«C(?.9)  ,WN0Z-C(30) 
C(21)=((104.+(9.65*EXPRN))*(ATA)**1.5)*PAMX/1000000. 

C  ( 2  ? )  =»  773000 .  /  SDMST 
C(?3)«C(.?1)*C(2?) 

C(?4)=1.53*ATA**.9 
C(?5)=»1?8800  ,/SDNTH 
C(?6)=C(74)*C(?5) 

C(?7)=.00?19*ATA*(EXPRN-4.)*TBMX*(1./CD)**1.7*PAKI!**.  9/1000000. 
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.11 749 
jnov 
JPM) 
JIW 
JIW 
JI9^ 
J?10'v 
,l???4 
J7774 
J??40 

J?33? 
J7  349 
J2452 
J757? 
J2609 
J?409 
J?409 
J?'7'> 

J77?*' 

J?73? 

J7900 

J7934 

J7944 

J?G30 

J794G 

J?090 

J3099 

j?r^ 

J3194 

J3204 

J?22a 

.13309 

J330? 

J344’1- 

J35'‘9 

J3632 

J37l$ 

J3799 

J3B72 

J395* 

.I'll  ?/•. 

j/iinu 
J47.49 
J<‘3?3 
JV::  2 
J4494 
J45W 
J4452 
J4736 
J4796 
J4976 
J'-vqRO 
J504? 
J5104 


c(?9)«4oooo./sdnxn 

c(?9)»c(?7)*c(?9) 

C(30)-C(23)+C(26)+C(?t») 

C  COMPONENT  LENGTH  CALCULATION 

FWDLT!“?.*RI>CHF*(  (EL0DF*ERATF+1  .  )/ERATF) 

AFTLT«?  .*RDCHA*(  (EL0DA*ERATA+1 .  t/ERATA ) 
XHOLT«((ATA/3.l4.159)**.q*((EXPRH)**i5-1,))/.3153 

C  MOTOR  WEIGHT,  MASS  FRACTION,  AMD  DELTA  VELOCITY  CALCULATION 
WCHR^V/CASE+VMMSL 

W I NRT  « ( WCHB+C ( 70 )  +VI I  GT+C ( 30  )  )*1 .0? 

WPROP**V/A+WF 

XMFRC®  ( WA+WF )  /  ( V/PR OP+W I  NRT  ) 

DE L V* ( 3 ? . 1 74 ) * ( X I  PR / WPR  OP ) *L OG  ( ( W I NRT+WPR OP ) / ( W I NRT ) ) 

WMTR»W !  NR  T  +V7PR  OP 

C  IF  GRAIN  CONFIG  IS  TANDEM  ADD  ALL  LENGTHS.  IF  GRAIN  CONFIG  IS 

C  PARALLEL  CHECK  FOR  LONGER  OF  FWD  AND  AFT  CKAM.fcLRS  AND  ADD  NOZZLE 

IF(MCAG-3)300,301 ,301 
3C0  RLGT»FV/!)l.T+AFTLT+XNOLT 


GO  TO  304 

30 1  I  F(  FV/DLT-A  FTLT  ) 30?  ,  30?. ,  303 
30?  RLGT«AFTLT+XNOLT 
GO  TO  304 

303  RI.GT«FVjDLT+XNOI.T 

304  I  FfRDCHF-RDCHA  'I3QJ;,30G  ,304 
30G  RDIA=? .*RDCHA+?.*TCYl>? ,*TINS 

GO  TO  307 

304  RD|A«? . * { RDCH F +FT C YL*FT  INS) 

307  ELODM^RLGVRDIA 

I  F( SENSE  SWITCH  1)300,309 

309  I F( SENSE  SV/ITCH  ?)  ?.  11,310 

311  PRINT  511, C(1), CONI ,  PSTAR,  FV/DLT,  A  FTLT  XMOLT 

310  PR  I HT50 1 , SDN ST , SDNTH , SDNXN , SDCHA , SDCHF, X i PR 
PRINT  50? , TE  MA , TE  MF , PA  MX, P  FMX , PA  MM 

PR  I  MT503  EI.ODA ,  ELODF  ,ERATA ,  ERATF  ,D  I A  A  ,D  I A  F 
PR  I  NT  qo  4  ,  CH  nil'.l,  PRC  ON ,  V/A ,  V/F ,  T  RMX ,  E  XPR  N 
PRINT  qnq,ATA,ATF,ASA,ASF,RO!NA  JROINF 
PRINT  q06,FLTP| ,WMAG,WMFG,ALPA,ALPF 
PRINT  q07,SFIG,V7M1G,EMG,GSIG,C0MIG,XMFIG 
PRINT  509,ENPA,ENPF,CPA,CPF,THC0A,THC0F 
PR  I  NT  509 , W I GT ,W I GP , W I G ,RH0A ,RHOF, CD 
PRINT  510,V/IHSL,V/IIIFC,WINAC,WINCY,Wi:!3,itCFG 
PRINT  5 1 1 , F I  NFC , F I NAC , F I NC Y , FW I  NS 
PR  I  NT  q  1  ?.  ,V/CASE  ,WFCL  ,WACL  ,WCYL ,  V/CHBT ,  MCAG 
PR  I  NT  513,  FV/FCL ,  FWACL ,  FWCYL  ,  FCH BT 
PR  I  NT  q 1 4,TCYL , T I  NS , FTC YL , FT  I  MS .ROCHA .ROCHF 
PRINT  515,C(?0),C(9).C(10).C(13),C(11)fC(1?) 
PRINT  516,C(19),C(18),C(17),C(16),C(15),C(14) 
PRINT  517,C(30),C(?.3),C(24),C(?.5)),EXPRN 
309  PRINT  5in,WMTR,RDIA,RLGT,XMFRC,DELV,EL0DM 


PUNCH  700 , XMFRC ,DELV,RLGT ,RDI A 
70  FORMAT  (7.9HS INGLE  CHAMBER  MOTOR  FOLLOWS) 
700  FORMA T(30H 
501  FORMAT (4H  100,4F11.0) 

507  FORMAT  ( 4H  700,6F11.4) 

503  F0RMAT(4H  300 ,6F1 1.5) 


,F11.3,ni.O,?F11.2) 
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J5195  SO''.  FORMAT  (Ml  *i00 ,  <>F1 1 .  ^ ) 

JC'TT  305  F0RMAT(M!  500, 5F  1 !  ,'n 

**?r«»o  so  a  format?  mi  aoo.sfu.s^ 

JS?r»r  S07  FORMAT'  Ml  700,6FH.M 

jsni.'i  sm  format  (Mi  noo,6Fii.*i\ 

JS'-?^  SOO  FORMAT  ('ll!  MOO/iFII.Sj 

>5531  Sin  FORMATf  M!1000,5F1 1.3.111) 
J5390  S11  F0RMA.T( '»H  1100 /iF1 1.31 

!5$G?  Sir  FORMATf  MM  ?00,5F1 1 .3.  Ill ) 

J S 7? >’■  S 1 3  FORMAT ( Ml  1 300 , 6f  1 1 . 3 ) 

.15733  fi'v  F0RMAT(nnroo,f>Fn.'O 

J53--0  515  FORMAT  (Mil  500  yfiF11.M 

J5910  m3  FORMAT  (Mil  600, F, F1 1  JO 

J507?  917  FORMATf MI1700.6 FI  1.M 

1 9.0 ^ ‘1.  511  FORMATf  Mil  100  ,f>F11 .4) 

,13035  OOf'TROL  10? 

.13103  GO  TO  10 

J6l 1 9  EMM 

PROOF?? I  MG  COMPLETE 
<1ART 
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APPENDIX  B  -  PRESENTATION  OF  DATA 


A  detailed  presentation  of  data  from  this  study  is  given  herein  in 
Figures  B-l  through  B-66.  For  ease  of  reference,  the  figures  are 
presented  by  'phase*  and  are  arranged  in  sequence.  Figures  B-l  through 
B-30  present  the  results  from  Phase  I;  the  remaining  .figures  are  for 
Phase  II. 

The  plots  for  Phase  I  are  presented  in  the  following  sequence: 

•  Mass  fraction  versus  minimum  thrust  at  an  I  ■  of  265 

lb  -sec/lb  for  1,  10,  and  20  starts  (Figures^B- 1,  B-2, 
and  B-3,  respectively^;  at  an  I.  of  280  lb^-sec/lb^  for 
1,  10,  and  20  starts  (Figures  B-4,  B-5,  and  B-6,m 
respectively);  and  at  an  I.  of  300  lb,-sec/lb  for  1,  10, 
and  20  starts  (Figures  B-T,  B-8,  ana  B-9).  T&'or  each 
plot,  curves  are  shown  for  total  impulse  values  of  104, 

10s,  and  106  lb  -sec  at  constant  thrust  ratios  of  1,  5,  and 

20.  f 

•  Mass  fraction  versus  total  impulse  for  burn  times  of  20, 

50,  200,  and  500  sec  (Figures  B-10  through  B-13, 
respectively).  All  curves  are  for  10  starts  and  constant 
specific  impulse  of  280  lb^-sec/lb^. 

•  Delta  velocity  versus  minimum  thrust  at  an  I.  of  265 
lb  -sec/lb-  for  1,  10,  and  20  starts  (Figures  B- 14, 

B-15,  and^-lb,  respectively);  at  an  I  of  280  lb^-sec/ 
lb  for  1,  10,  and  20  starts  (Figures  B-17,  B-18,  and 
B^J9,  respectively); and  at  an  I.  of  300  lb,-sec/lb  for 
1,  10,  and  20  starts  (Figures  B-20,  B-21,  and  B- 22, 
respectively).  For<  each  plot,  curves  are  shown  for 
total  impulse  values  of  10  ,  ldB,  and  108lb^-sec  at 
constant  thrust  ratios  of  1,  5,  and  20. 

•  Delta  velocity  versus  total  impulse  for  burn  times  of 
20,  50,  200,  and  500  sec  (Figures  B-23  through  B-26, 
respectively).  All  curves  are  for  10  starts  and  constant 
specific  impulse  of  280  Ib^-sec/lb^. 
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Length  versus  minimum  thrust  for  1^  of  265  and  300  lb  -sec/lb 
(Figures  B-27  and  B-28,  respectively;.  Curves  of  constant  m 
configuration  are  shown  on  each  plot. 

Diameter  versus  minimum  thrust  for  I^p  of  265  and  300  lb  -sec/lb 
(Figures  B-29  and  B-30,  respectively).  Curves  of  constant  m 

configuration  and  constant  total  impulse  values  (104,  106 , 


5  X  10B,  and  10®  lbf-  sec)  are  shown  on  each  plot. 
The  plots  for  Phase  II  are  arranged  as  follows: 
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Mass  fraction  versus  minimum  thrust  at  thrust  ratios  of  1,  5, 
and  20  (Figures  B--3J.,  B-32,  and  B-33,  respectively).  Curves 
for  three  materials  are  shown  on  each  plot. 

Delta  velocity  versus  minimum  thrust  at  thrust  ratios  of  1,  5, 
and  20  (Figures  B-34,  B-35,  and  B-36,  respectively).  Curves 
for  three  materials  are  shown  on  each  plot. 

Mass  fraction  versus  minimum  thrust  for  aft  propellant  burning- 
rate  constants  of  0.  0002,  0.  0004,  and  0.  0010  (Figures  B-37, 
B-38,  and  B-39,  respectively).  Curves  for  constant  forward 
propellant  burning  rates  of  0.  002,  0.  0008,  and  0.  0032  are 
shown  on  each  plot. 

Delta  velocity  versus  minimum  thrust  for  aft  propellant 
burning-rate  constants  of  0.  0002,  0.  0004,  and  0.  0010  (Figures 
B-40,  B-41,  and  B-42,  respectively).  Curves  for  constant  for¬ 
ward  propellant  burning  rates  of  0.  0002,  0.  0008,  and  0.  0032 
are  shown  on  each  plot. 

Length  versus  aft  propellant  burning-rate  constant  at  minimum 
thrust  values  of  200,  500,  1000,  and  5000  lb^  (Figures  B-43 
through  B-46,  respectively).  Curves  for  constant  forward 
propellant  burning  rates  of  0.0002,  0.0008,  and  0.0032  are 
shown  on  each  plot. 

Diameter  versus  aft  propellant  burning-rate  constant  at 
minimum  thrust  values  of  200,  500,  1000,  and  5000  lb^  (Figures 
B-47,  through  B-50,  respectively).  Curves  for  constant 
forward  propellant  burning  rates  of  0.  0002,  0.  0008,  and 
0.  0032  are  shown  on  each  plot. 
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•  Mass  fraction  versus  thrust  ratio  at  forward  propellant  burning- 

rate  exponents  of  0.  6  0.  8,  and  0.9  (Figures  B-51,  B-52,  and 

B-53,  respectively).  Curves  for  constant  aft  propellant  burning 
rate  exponents  of  0.  8,  1.  0,  and  1.  1  are  shown  on  each  plot. 

•  Delta  velocity  versus  thrust  ratio  at  forward  propellant  burning- 
rate  exponents  of  0.  6,  0.  8,  and  0.  9  (Figures  B-54,  B-55,  and 
B-56,  respectively).  Curves  for  constant  aft  propellant 
burning-rate  exponents  of  0.  8,  1.  0,  and  1.  1  are  shown  on  each 
plot. 

•  Length  versus  thrust  ratio  at  forward  propellant  burning-rate 
exponents  of  0.  6,  0.  8,  and  0.  9  (Figures  B-57,  B-58,  and 
B-59,  respectively).  Curves  for  constant  aft  propellant 
burning-rate  exponents  of  0.  8,  1.  0,  and  l.  l  are  shown  on  each 
plot. 

•  Diameter  versus  thrust  ratio  at  forward  propellant  burning- 
rate  exponents  of  0.6,  0.8,  and  0.  )  (Figures  B-60,  B-61,  and 
B-62,  respectively).  Curves  for  constant  aft  propellant  burning- 
rate  exponents  of  0.  8,  1.  0,  and  1.  1  are  shown  on  each  plot. 

•  Mass  fraction  versus  minimum  thrust  (Figure  B-63).  Curves 
for  constant  theta  values  of  2,  3,  and  4  are  shown. 

•  Delta  velocity  versus  minimum  thrust  (Figure  B-64).  Curves 
for  constant  theta  values  of  2,  3,  and  4  are  shown. 

•  Length  versus  minimum  thrust  (Figure  B-65).  Curves  for 
constant  theta  values  of  2,  3,  and  4  are  shown. 

•  Diameter  versus  minimum  thrust  (Figure  B-66).  Curves 
for  constant  theta  values  of  2,  3,  and  4  are  shown. 
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Figure  E  1  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  265  Lb.-Sec/Lb  and  1  Start 
sp  f  m 
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Figure  B-3  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  265  Lb,-Sec/Lb  and  20  Starts 
sp  f  m 
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Figure  B-4  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  280  Lb.-Sec/Lb  and  1  Start 
sp  f  m 
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Figure  B-5  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  280  Lb,-Sec/L.b  and  10  Starts 
sp  f  m 
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Figure  B-6  -  Mass  Fraction  versus  Minimum  Thrust  for 
“  I  =  280  Lb  -Sec /Lb  and  20  Starts 
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ass  Fraction  versus  Minimum  Thrust  for 

I  a  300  Lb.-Sec/Lb  and  1  Start 
sp  z  m 
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Figure  B-8  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  300  Lb,-Sec/Lb  and  10  Starts 
sp  f  m 

CONFIDENTIAL 


B-ll 


Minimum  Thrust  (lb  ) 


CONFIDENTIAL 


AFRPL-TR-65-209,  Vol  II 


uowtxg,  ■•<*W 


B- 12 


Figure  B-9  -  Mass  Fraction  versus  Minimum  Thrust  for 

I  =  300  Lb,-Sec/Lb  and  20  Starts 
sp  f  m 
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Figure  B-10  -  Mass  Fraction  ve-sus  Total  Impulse  for 
-20-Sec  Burn  Time  . 
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Figure  B-ll  -  Mass  Fraction  versus  Total  Impulse  for 
50-Sec  Burn  Time 
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Figure  B-12  -  Mass  Fraction  versus  'I'otal  Impulse  for 
200-Sec  Burn  Time 
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Figure  B-13  -  Mass  Fraction  versus  Total  Impulse  for 
500-Sec  Burn  Time 
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Figure  B- 14 


Delta  Velocity  versus  Minimum  Thrust  for 

I  =  265  Lb  -Sec/L.b  and  1  Start 
sp  f  m 
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Figure  B-15  -  Delta  Velocity  versus  Minimum  Thrust  for 

I  =  265  Lb,-Sec/Lb  and  10  Starts 
sp  f  m 
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Figure  B- 17 


B-  20 


Delta  Velocity  versus  Minimum  Thrust  for 
I  =  280  Lb,-Sec/Lb  and  1  Start 
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Figure  B-19  -  Delta  Velocity  versus  Minimum  Thrust  for 
I  =  280  Lb,-Sec/Lb  and  20  Starts 
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Figure  B-20  -  Delta  Velocity  versus  Minimum  Thrust  for 
I  =  300  Lb.-Sec/Lb  and  1  Start 
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Figure  B-21  -  Delta  Velocity  versus  Minimum  Thrust  fo 

I  =  300  Lb,-Sec/Lb  and  10  Starts 
sp  i  m 
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Figure  B-  22  -  Delta  Velocity  versus  Minimum  Thrust  for 

I  =  300  Lb.-Sec/Lb  and  20  Starts 
so  £  m 
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Figure  B-23  -  Delta  Velocity  versus  Total  Impulse 
for  20-Sec  Burn  Time 
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Figure  B-24  -  Delta  Velocity  versus  Total  Impulse 
for  50-Sec  Burn  Time 
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Figure  R-25  -  Delta  Velocity  versus  Total  Impulse 
for  200-Sec  Burn  Time 
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Fi  »’  re  B-26  -  Delta  Velocity  versus  Total  Impulse 
for  500-Sec  Burn  Trine 
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Figure  B-27  -  Length  versus  Minimum  Thrust  for 

I  »  265  Lb.-Sec/Lb 
sp  f  m 
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Figure  B-29 


Diameter  versus  Minimum  Thurst 


for  I  =  265  Lb  .-Sec /Lb 
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Figure  B-30  -  Diameter  versus  Minimum  Thrust 
for  I  =  300  Lb,-Sec/Lb 
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Figure  B-36  -  Delta  Velocity  versus  Minimum 


Thrust  at  Thrust  Ratio  of  20 
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Figure  B-37  -  Mass  Fraction  versus  Minimum  Thrust  for  Aft 
Propellant  Burning- Rate  Constant  of  0.  000  2 
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/lass  Fraction  versus  Minimum  Thrust  for  Aft 
^ropellant  Burning- Rate  Constant  of  0.  0004 
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Figure  B-39  -  Mass  Fraction  versus  Minimum  Thrust  for  Aft 
Propellant  Burning-Rate  Constant  of  0.  0010 
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Propellant  Burning-Rate  Constant  of  0.  0002 
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Figure  B-41  -  Delta  Velocity  versus  Minimum  Thrust  for  Aft 
Propellant  Burning-Rate  Constant  of  0.  0004 
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Figure  B-42  -  Delta  Velocity  versus. Minimum  Thrust  for  Aft 
Propellant  Burning-Rate  Constant  of  0.  0010 
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Figure  B-44  -  Length  versus  Aft  Propellant  Burning-Rate 

Constant  for  F  =  500  Lb, 
mxn  f 
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Figure  B-48  -  Diameter  versus  Aft  Propellant  Burning- 

Rate  Constant  for  F  =  500  Lb, 

mm  f  B-51 
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Figure  B-49  -  Diameter  versus  Aft  Propellant  Burning 
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Rate  Constant  for  F  =  1000  Lb, 
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Figure  B-51  -  Mass  Fraction  versus  Thrust  Ratio  for  Forward 
^  Propellant  Burning-Rate  Exponent  of  0.  6 
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Figure  B-52  -  Mass  Fraction  versus  Thrust  Ratio  for  Forward 
Propellant  Burning-Rate  Exponent  of  0.  8 
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Figure  B-53  -  Mass  Fraction  versus  Thrust  Ratio  for  Forward 
r  Propellant  Burning-Rate  Exponent  of  0.  9 
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Figure  B-55  -  Delta  Velocity  versus  Thrust  Ratio  for  Forward 


Propellant  .Burning -Rate  Exponent  of  0.  8 
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Figure  B-56  -  Delta  Velocity  versus  Thrust  Ratio  for  Forwar< 
Propellant  Burning-Rate  Exponent  of  0.  9 
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Figure  B-57  -  Length  versus  Thrust  Ratio  for  Forward  Propellant 
Burning -Rate  Exponent  of  0.  6 
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B-62 


Length  versus  Thrust  Ratio  for  Forward  Propellant 
Burning -Rate  Exponent  of  0.  9 
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Figure  B-61  -  Diameter  versus  Thrust  Ratio  for  Forward 
Propellant  Burning-Rate  Exponent  of  0.  8 

CONFIDENTIAL 


CONFIDENTIAL 


AFRPL-TR- 65-209,  Vol  II 


'igure  B-62  -  Diameter  versus  Thrust  Ratio  for  Forward 
Propellant  Burning -Rate  Exponent  of  0.  9 
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Figure  B-63  -  Mass  Fraction  versus  Minimum  Thrust  for  Theta 


Values  of  2,  3,  and  4 
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Figure  B-64  -  Delta  Velocity  versus  Minimum  Thrust  for  Theta 
Values  of  2,  3,  and  4 
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Figure  B-65  -  Length  versus  Minimum  Thrust  for  Theta  Values  of 
2,  3,  and  4 

CONFIDENTIAL 


B-68 


Minimum  Thrust  (lb  ) 


CONFIDENTIAL 

AFRPL-TR-65-209,  Vol  II 


N 


-O  «VJ  00  ©  vO  M 

to  .pi  p-J  <m  r>j  -*  — « 


(  'Uf )  aO^OOTBTQ 


igure  B- 66  -  Diameter  versus  Minimum  Thrust  tor  Theta  Values 
of  2,  3,  and  4 
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